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SUMMARY. 

Sound  is  basically  the  vibration  of  continuous  elastic  substances.  The  vibra¬ 
tion  takes  the  form  of  an  acoustic  wave  which  can  be  mathematically  analyzed  in 
the  same  manner  as  an  electric  wave  can  be  analyzed.  Sound  waves  generated  in  an 
acoustic  medium  and  confined  by  the  rigid  boundaries  of  a  tube  may  be  described 
in  much  the  same  manner  as  one-dimensional  transverse  waves  in  a  vibrating  string 
of  finite  length  or  by  two-dimensional  waves  in  a  membrane  fixed  along  its  edge. 

The  first  part  of  this  report  concerns  itself  with  the  development  of  the 
wave  equation  in  its  general  form.  For  this  purpose,  there  is  a  discussion  of 
the  Equation  of  State,  the  Newtonian  Mechanical  Equations  as  applied  to  an  ideal 
fluid,  and  the  Law  of  Conservation  of  Mass.  These  relationships  lead  to  a  second 
order  partial  differential  equation  known  as  the  three-dimensional  wave  equation. 
It  essentially  describes  a  wave  being  propagated  in  a  deformable  medium  in  three 
mutually  perpendicular  directions.  In  order  to  apply  the  wave  equation  to  the 
analysis  of  sound  waves  in  water  filled  tubes,  the  equation  is  converted  from 
Cartesian  to  cylindrical  coordinates. 

It  is  then  shewn  that  the  wave  equation  is  separable  and-that  its  solution 
may  be  represented  as  the  product  of  four  distinct  functions,  each  related  to 
only  one  variable.  Each  of  these  product  functions  is  solved  by  successively 
reducing  the  partial  differential  equation  to  several  ordinary  differential 
equations.  The  general  solution  to  the  wave  equation  is  then  shown  to  be  a 
complex  waveform  being  simultaneously  propagated  along  the  longitudinal  and 
radial  axes  of  the  tube  as  well  as  in  the  ©-space.  The  net  result  is  a  wave 
which  is  propagated  with  a  screw-like  motion  along  the  longitudinal  axis  of 
the  tube.  The  wave  may  be  either  progressive  or  stationary,  depending  upon  given 
boundary  conditions.  Assuming  that  the  walls  of  the  tube  in  which  the  sound  wave 
is  being  propagated  form  a  perfect  reflector  and  that  the  waves  must  be  finite 
at  the  origin  of  the  tube,  it  is  shown  that  the  wave  equation  reduces  to  a  plane 
wave,  that  is,  one  which  is  propagated  along  the  longitudinal  axis  of  the  tube 
only.  It  is  further  shown  that,  as  the  radial  waves  form  a  discontinuous 
series  of  frequencies,  a  certain  minimum  frequency  must  not  be  exceeded  in  order 
to  maintain  plane  monochromatic  waves. 

The  final  part  of  this  report  deals  with  the  application  of  theory  to  sound 
waves  in  .water  filled. tubes.  The  acoustic  impedance  of  a  material  placed  at 
one  end  of  the  tube  is  shownto  be  a  complex  quantity  which  is  a  function  of 
known  constants  of  the  water,  sound  wave  and  tube.  It  is  also  a  function  of  the 
phase  angle  between  incident  and  refleqted  rays  at  the  interface  of  the  water 
and  material  and  of  the  amplitude  ratio  of  reflected  to  incident  pressure  waves 
(pulse  tube  method);  or  the  ratio  of  minimum  to  maximum  pressure  and  the  dis¬ 
tance  cf  the  pressure  nodes  from  the  interface  (standing  wave  method).  Pressure 
ratios  and  phase  angles  are  measureable  quantities.  Once  determined,  the 
acoustic  impedance  may  be  calculated  or  graphically  determined  from  Smith  Charts. 
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1 «  Definition  of  k  Fluid.  Variables  Describing  Continuous  Deformations-1- f 2 


Acoustic  waves  are  generally  generated  byieontanudLsand  oscillatory  changes  in 
the  physical  properties  of  matter  (pressure,  density,  etc.).  They  give  rise  to 
audible  sounds  when  their  frequencies  vary  from  about  20  cps  to  15,000  cps. 

There  are  also  ultrasonic  frequencies  above  15,000  cps  and  subsonic  frequencies 
below  20  cps.  The  manner  in  which  a  material  transmits,  absorbs  or  reflects 
acoustic  waves  is  closely  connected  with  its  internal  structure.  Yet  the 
relationships  between  the  elemental  structure  of  materials  and  the  values  of 
their  empirical  constants  has  not  been  calculated  to  date. 

The  way  continuous  deformations  give  rise  to  acoustic  waves  is  a  consequence 
of  the  Newtonian  mechanical  equations  for  fluids,  plus  certain  supplementary 
conditions  and  approximations . 

Let  us  assume  that  we  are  dealing  with  a  quasi -continuous  medium  like  air, 
water  or  steel.  The  first  two  media  are  called  "fluids"  because  of  their 
propensity  to  flow.  But  the  theory  of  relativity  taught  us  that  between  the 
three  substances  there  is  only  a  difference  of  degree,  a  solid  under  extreme 
conditions  of  pressure  and  temperature  being  allowed  to  creep  like  a  liquid. 

Thus,  any  substance  can  be  called  a  "fluid".  Fluid  properties  are  defined  by 
pressure  F,  mass  density  a,  specific  weighty,  viscosity  u,  surface  tension  6  , 
and  modulus  of  elasticity E  .  For  a  low  Mach  number  M  e  ~  (ratio  of  velocity 
of  a  fluid  v  with  respect  to  a  given  coordinate  system  to  sound  velocity  c  in 
the  fluid),  both  hydro -dynamics  and  aerodynamics  may  be  treated  in  the  same 
manner.  Compressibility  becomes  important  and  must  be  taken  into  account  only 
when  the  Mhch  number  is  over  about  «. 

In  a  fluid  we  can  define  at  each  point  P  at  time  t  a  pressure  P,  a  density  j 
and  a  velocity  with  components  u,  v,  w  with  respect  to  a  rectangular  coordinate 
system,  Which  is  the  local  velocity  of  a  particle  of  the  fluid  passing  a  point 
P  (ijS,  y,  i)  at  time  t  and  which  could  be  very  different  from  the  average  velocity 
of  the  whole  fluid.  The  five  quantities  P,  o,  u,  v  and  w  are  in  general  functions 
of  the  coordinates  X,  y,  %  and  time  t.  In  order  to  determine  the  behavior  of  each 
element  or  particle  of  the  fluid  when  t  varies,  we  must  know  the  values  of  these 
five  quantities  at  any  time.  Therefore,  we  must  have  five  relations  between 
these  five  quantities.  The  usual  way  to  get  them  is  to  write  five  differential 
equations  determining  them  in  terms  of  five  arbitrary  initial  constants  P0#  a  0> 
u0J  v0  and  w0  representing  the  values  of  f,  u,  v  and  w  at  time  t  =  t0.  ■ 
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1  2 

2.  Equation  of  State  ’ 

The  first  relation  we  need  may  be  derived  from  the  equation  of  state*  In 
any  fluid  there  is  a  relationship  between  pressure  P,  density^  and  absolute 
temperature  T : 

1)  y,  (P,o,  T)  -  0 

For  an  ideal  gas  at  constant  temperature  (1)  takes  the  form 

2)  Py  =  RT 

> 

where  F  is  the  absolute  pressure,  v  is  the  specific  volume,  T  the  absolute 
temperature,  and  R  the  universal  gas  constant.  On  setting  v  =  l/a  we  get 
P  =oRTj^  In  the  adiabatic  case  equation  (?)  is  non-linear  and  takes  the 
form  Pv°*=  constant  or 


3)  P  ■  g  (a)  or  a  =  f  (P) 
whereof  is  the  gas  constant. 

For  perturbations  of  small  amplitudes  (which  is  usually  the  case  for  acoustic 
waves  if  we  disregard  shock  waves)  equation  (3)  can  be  placed  by  a  linear  rela¬ 
tion,  obtained  by  retaining  the  first  two  terms  oC  a  ^aylor  expansion  of  (3)  and 
disregarding  higher  order  terms. 


so  that  A  P  =  P-Po  = 
have  from  equation  (li)'d 


Ao 


Since  the  velocity  of  sound  c  we 


5)  P  =  Po  +  c^  Ao  '  . 

This  equation  may  be  regarded  as  the  equation  of  state" for  perturbations  of 
small  amplitude  in  ideal  gases.  The  equation  of  state  fo,r  real'  gases  at  constant 
temperature  would  closely  approximate  Van  der  Waal's  equation 

6)  . (P  +  a/V2)  (v-b)  =  RT 


where  P  and  V  are  the  pressure  and  volume  and  a  and  b  are  molecular  constants. 
Real  gas  equations  of  state  crtn  also  represent  the  behavior  of  liquids  within 
a  limited  range  of  pressures  and  temperatures.  But,  for  a  liquid  in  bulk,  the 
equation  of  state  does  not  contain  F  as  any  liquid  is  very  nearly  incompressible. 
Density  a  depends  only  on  temperature  T  and  we  have  between  certain  limits 
0  eo0  (1  +  KT)"1,  K  being  the  coefficient  of  dilation  of  the  liquid  due  to  heat 
and  o0  its  density  at  T  =  0.  For  most  cases,  water  may  be  regarded  as  incom- 
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pressible,  but  not  as  far  as  sound  waves  are  concerned.  In  other  words,  water 
appears  incompressible  for  static  pressures  (between  certain  limits).  The 
passage  of  a  sound  wave  through  water  is  really  a  pressure  wave  with  a  small 
variation  of  density  approximately .expressed  by  Ap  =  lq  ,  the  large  value  of 
ir  associating  a  small  variation  of  density  to  to  a  large  variation  of 
pressure  Ap. 
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3.  Ftefinition  of  External  Forces 

The  four  remaining  equations  are  derived  from  the  three  equations  of  motion 
(Newtonian  mechanical  equations  applied  to  fluids)  and  the  law  of  mass  conserva¬ 
tion. 


Suppose  that  at  time  t  a  fluid  F  occupies  a  volume  V  bounded  by  a  surface  S, _ 

satisfying  the  usual  conditions  of  smoothness  and  uniqueness.  The  external  forces 
acting  on  the  elements  of  this  fluid  are  of  two  kinds:  those  acting  on  the 
volume  elements  occupied  by  F  and  those  acting  on  the  surface  elements  immersed 
in  or  limiting  F. 

a.  The  external  forces  in  a  volume,  like  gravity,  electric  field,  hydrostatic 
pressure,  etc.,  .  .  .  are  defined  in  terms  of  the  element  of  mass  dm  of  the  fluid 
F  contained  in  the  volume  dV,  so  that  dm  »  odV,  o  being  the  density  of  the  fluid  F 
at  a  point  P  in  the  volumetric  element  uV.  These  elementary  forces  fv  are  con¬ 
sidered  to  be  proportional  to  dm  and  therefore  are  defined  by  the  equation 

7)  dfv  -  Kydm  *»  K^dV 

Kv  being  a  proportionality  constant  equal  to  the  force  per  unit  mass. 

b.  There  are  also  forces  acting  through  areas  immersed  in  or  limiting  F, 
such  as  pressure  exerted  by  the  walls  of  the  container  confining  F.  Let*  dS  be 
an  elemental  area.  The  surface  forces  fs  are  considered  to  be  proportional  to 
the  area  upon  which -they  act  and  thus  we  may  write 

8)  dfa  -  K3  dS 

where  Ks  is  a  proportionality  constant  equal  to  the  force  per  unit  area. 
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4.  Mechanical  Equations' 


We  may  now  write  the  mechanical  equations  of  motion  for  a  fluid  F  corresponding 
to  the  Newtonian  equations  for  any  material  system.  This  is  simply  an  applica¬ 
tion  of  Newtonian  mechanics  to  continuous  media.  Consider  a  point  P  in  the  fluid 
with  mass  m,  velocity  (u,  v,  w)  submitted  to  a  force  (Fx,  F  ,  F  )  with  respect  to 
the  (x,  y,  z  )  coordinate  system.  The  Newtonian  equations  of  motion  are: 


9)  _  du 

F*=mZt 


r-  i 

V  *"31 


>  F*  =  ^  3BE? 


If  X,  Y,  Z  are  components  of  the  force  per  unit  mass  acting  on  the  fluid  F 
at  P  then  - 

10>  X*  <fV=  YsdLV= 

are  components  of  the  total  force  acting  on  the  mass  dm.  If  J  is  the  total 
acceleration  of  a  particle  with  coordinates  x,  y,  z  and  mass  dm,  then  tne  com¬ 
ponents  of  J  are 

T  T  ,  T  cL2Z 

v  dltY  ?  djt* .  'l^*^^The  elementary  inertial  force  dm  J  applied 

to  the  particle  has  the  components  dm  Jx,  dm  J  and  dm  Jz»  Integration  of  (10) 
yields  the  components  of  the  total  force  acting  on  the  volume  V  occupied  by 
fluid  F  at  time  t. 

■»  SB^ . 

Similarly  the  total  acceleration  for  the  same  volume  V  at  time  t  has  the  components 

JJl  J^V 

giving  rise  to  inertial  forces  fff  TyAcLV the 


giving  rise  to  inertial  forces  Jjf^dZdV^  fjf  Jy^ctV  ^e  the 

fluid  is  in  a  State  of  equilibrium,  these  volume  forces  must  be  compensated  by 
the  total  force  through  the  surface  S  limiting  the  volume  V.  This  surface  force 
has  the  components 


Ji^dS  ,  ffs  T *  ds,  ;;s t2  ds 

the  mechanical  equations  of  motion  become 

JJL  ^(X-J ;)AV-jT  Ty.d  S  =  O  _ 


u,  JUv  A(y-Ly) dv-ffs  Ts/d$  =  c> 
fjfv  *  Cz- J2)  dv-  SL  ~rz  cLs  =  o 

We  may  also  say  that  the  sum  of  all  external  forces  must  be  equal  to  the 
inertial  force. 
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Euler’s  Equations  for  an  Ideal  Fluid 


Equations  (l)i)  become  simpler  when  we  consider  the  case  of  an  ideal  (or 
perfect)  fluid  for  which  the  viscosity  is  negligible  and  the  shear  modulus  is 
zero.  Usual  fluids  like  water  as  far  as  acoustic  waves  are  concerned  aj  proach 
ideality  sufficiently  to  make  use  of  this  simplification.  Jt  is  in  agreement 
with  Prandtl’s  hypothesis  according  to  which,  with  fluids  of  low  viscosity, 
the  effects  of  viscosity  are  limited  to  a  narrow  region  along  the  boundaries. 
In  many  real  cases,  fluid  friction  is  so  small  that  the  hypothesis  of  ideality 
is  sufficiently  accurate.  For  an  ideal  fluid,  the  surface  forces  through  dS 
are  all  normal  to  dF  and  directed  in  the  positive  direction.  Denoting  asc(, 

(3,  Y  the  direction  cosines  of  the  outer  normal  to  S  at  dS,  we  get  fer  the 
components  of  the  forces  Tx,  T,r  and 

p, 

15) 


fV 


Pxct 


■X 


dP-w 


Pv  P 


T„  - 


■I 


dP, 


pzr  = 


1 


dP„ 


where  Px,  Py  and  P2  are  the  components  of  pressure  in  the  x,  y  and  z  directions. 
dpx  dPy  "(if)  dy’  dp*  "(»)dz’  we  set_/dp*  “/(a*)  *c  ■ 

~ _ _ x  n  __  f 


Since 


etc. 


f  •  ♦ 


consequently, 


16)  Tx  dF  =  JJ  rxc(dS  =  JJ  dPxj  dS  0x 

Substituting  in  equations  (lh)  we  obtain 


dS  ■ 


fffv  •)  dV"° 


These  equations  must  hold  for  iny  volume  whatsoever,  in  particular  for  an 
arbitrary  infinitesimal  volume  dV.  This  requires  that  the  integrand  vanish. 
Therefore 


18)  O  (X-Jx)  -  =  0 

o(Y"Jy)  -fy  =  0 

o(Z-Jfe)  -  &  =  0 
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Now,  let  us  follow  the  motion  of  a  certain  particle  of  the  fluid  in  terms 
of  time  t.  When  passing  the  point  P  at  time  t  it  has  the  velocity  (u,  v,  w) 
and  obviously  the  values  of  these  components  depend  on  coordinates  x,  y,  2  of 
the  particle.  In  turn,  the  latter  depend  on  time  t,  since  the  coordinates  of 
the  particle  vary  when  time  increases.  Thus,  u,v  and  w  are  functions  of  x,  y,  z 
and  c 


19)  dx/dt  =  u(x,  y*  z,  t) 

dy/dt  =  v(x,  y,  z,  t) 
dz/dt  =  w (x ,  y,  z,  t) 


and  x 
20) 


=  x(t),  y  =  y(t),  z  = 
T 

“  dli?*  dc t  ~ 

a?y  dtv  _ 

Jy =  dt2  d± 

r  _  dw 

"  dtx  "  dfc 


z(t).  Consequently  the  acceleration  becomes 


using  the  well  khown  formulas  of  partial  derivatives  and  the  chain  rule. 


Replacing  in  (18)  Jx,  Jy; 


J_  by  their  values  (20) ,  we  obtain  the  Euler  equations 


21) 


JL 

_L 

A 


dP 


a p 

dz. 


*Z- 


6t 


-  OL 


'Which,  when  integrated  give  three  conditions  to  be  satisfied  by  the  variables 
u,  v,  w.  Together  with  the  equation  of  state  (1)  taken  in  the  adiabatic  form 
not  involving  temperature  T,  i-_e;;  pv°^  =  constant,  since  temperature  has  no 
time  to  vary  during  the  processes  we  are  dealing  with,  we  obtain  four  conditions. 
We  still  need  another  equation  to  determine  the  variables  ? t  & ,  u,v  and  w.  It 
is  given  by  the  Law  of  Conservation  of  Mass;  and  the  Continuity  Equation, 
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6.  Continuity  Equation 

Consider  a  particle  of  the  perfect  fluid  (ideal  frictionless  medium)  with  mass 
dm  occupying  the  volume  dV  at  time  t  and  having  the  coordinates  x,  y,  z.  We  have 

dm  =  odV  where  0  is  tjjie  densitv  at  the  point  P(x,  y,  z)  at  time  t.  When  t  becomes 

t*  =  t  +  dt  is  at  Q(x  ,  y  ,  z  )  with  coordinates 

22)  x1  =  x  +  u  dt 
y'  =  y  +  v  dt 
a*  =  z  +  w  dt 

-/  1 

Th«  volume  occupied  by  the  same  particle  becomes  dv  ,  its  density  0  at 
time  t  but  the  mass  remains  the  same  so  that 

23)  dm  =  odV  =  <r'd  V* 

Since  a  is  a  function  of  x,  y,  z  and  t  we  may  write  a  “  f(x,  y,  z,  t)  so  that 

2lj)  or'  =  f(x*,  y' ,  z  ,  l  ) 

=  f(x  +  udt,  y  +  vdt ,  z  +  wdt,  t  +  dt) 

Using  a  Taylor  series  expansion  limited  to  the  second  term  we  have 

2^ )  cr#  =  a  +(\ ^ju  +(^)v  w  +  dt 
^3xj  [dyj  [dz/  3t/ 

On  the  other  hand  dV*  =  dx;  dy*  dz^ ,  The  term  dx*  dy*  dz^  may  be  computed 
by  differentiating  formulas  (22)  which  yield 

26)  dx*  =  dx  +  du  dt 

dyr  =  dy  +  dv  dt 

y 

dz  ®  dz  +  dw  dt 

Partial  differentiation  of  (22)  or  division  of  (26)  by  dx,  dy,  dz  gives 


t 

a*_ 

=  i+( 

a jl. 

dt 

$L.  J 

>\<it 

ax 

,dx) 

ay 

lay/ 

dz  \ 

ax 

11 

m. 

dt 

a_zl  = 

dy 

i+(in  dt 
\dy) 

$L_  J 
dz  1 

te)« 

J&l' 

jdx 

=/( 

^ t 

Wm.mm  «  1  +  J 

(l)dt 

ax 

\dxi 

oy  \i 

3z  i 

12 


Substituting  in  dV 
(dt)^  we  have 
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dx  dy  dz  and  neglecting  terms  containing  (dt)  and 


*v. 


From  (25)  we  deduce,  in  the  same  approximation,  using  (25)  and  (28) 


dV  [!  +  ( 


c)U,  dv  duJ 

6^  dy  dz 


)<&]  - 


cLV 


Neglecting  the  (dt)  term  (29)  becomes  . 

^av+-  At  av[4§+^  +  sl]+[(lf)^K^)v+(^)w+dt 

which  reduces  to 


so)  r^u.  <^V  dwl  . 

£S'U+a74-dTj  + 

(30)  can  be  written 


(MY*  (M) 


v  4- 


az)w  +  at 


315  ^ +  =o 

dt  dy.  T  ^7  « 

which  is  the  continuity  equation.  It  represents  the  conservation  of  mass  in  a 
differential  form.  Indeed,  integrating  (31)  over  volume  V  we  obtain 

By  Gauss’  theorem  for  a  vector  If  where 

~  -  (Kx,  Ky»  Kz)  we  have  JL  Kn  dS  =  f£fv div  *K  dV,  which  means  that  the  surface 
integral  of  the  normal  component  of  K  over  the  surface  S  .equals  the  volume 
integral  over  the  volume  V  bounded  by  S  where  div  K  =  ^Kz 

rrr/^m  Uv  tew\  ~st'  f  ^  +  ' 

The  integraljjj  ydV  °f  (32)  can  tJliere:fore  be  written  in  the 
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form  of  a  surface  integral,  surface  S  bounding  volume  V.  Denoting  as  K  the 
vector  with  components  au,0  v,  andow,  ffj  div  *K  dV  =  Jjf  Kn  dS.  If  S  is 
the  surface  bounding  the  fluid,  no  amount  of  fluid  flows  through  this  surface 
and  as  K  is  the  vector  through  this  surface  then  Kn  dS  must  be  zero,  where 
Kn  is  the  component  of  K  normal  to  the  surface.  The  integral 

III  (l°)dv  of  (3J) 

may  be  written  as 


U/// 


=  sl¬ 
at 


dm  =  |E. 

at 


'  v  s 

Since  the  second  integral  of  (3)  must  be  zero  then  (32)  reduces  to  =  0  and 


the  mass  is  conserved, 
inside  ?, 


This  implies  that  no  fluid  is  created  or  destroyed 
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7.  Force  Potential  and  Velocity  Potential 


The  Euler  equations  (21)  may  be  put  into  a  simpler  form  when  three  supple¬ 
mentary  restrictions  are  valid. 

First,  assume  that  the  external  volume  ?orce  X,  Y,  Z  derives  from  an  energy 
potential  V,  so  that 

33)  y  —  \/_  7  _  — 

5  /-  ay  1  L  <^z 

Second,  there  is  no  rotational  motion  inside  the  fluid  so  that  the  fluid 
velocity  must  derive  from  a  potential  $ 


w  =  - 


34)  V=-^  w=-^ 

The  negative  sign  indicates  that  an  increasing  energy  potential  V  causes  the 
velocity  of  the  fluid  to  decrease. 


Third,  the  equation  of  state  is  valid  in  the  adiabatic  formd~=  f  (P)  solved 
with  respect  to<£\ 

When  these  three  conditions  are  satisfied  we  may  define  a  new  function 

r  dip 


Gl’V-f 


Introducing  .his  equation  into  the  Euler  equations  (21)  and  replacing  X,  Y,  Z 
by  (33)  and  u,  v,  w  by  (34)  respectively,  we  obtain  for  the  x  component 

36)  =  _  X 

cKf  o>Y 

and  we  get  similar  formulas  for^^/^yand  s<%z  ,  replacing  x  in  (36)  by  y  and  z 
sucessively,  As  ^  is  a  uniform  function  of  x,  y,  z  and  t,  we  have 

crtl  axf  4-  f 

— »  _  — ±.  r  t — rr,  .  eTc  for  x ,  y,  z  „  t . 

o>Y^y“  c^Y^Sc  )  c)tc^Y  d\f.dt.) 

Thus,  the  right  hand  side  of  (36)  can  be  written  in  the  form 


c> 


M+-1 

oJ't  2. 


_Caf-)  +  (<37J+C3 
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consequently,  we  may  write  (36)  in  the  form 


The  quantity  within  the  outer  bracket, 


/  j  2  ^ 

Q  4-  &*  _  +  /&A 

at.  2  Max)  [dyj  \3z) 

must  therefore  be  independent  of  x  as  its  partial  derivative  with  respect  to  x 
is  zero.  Similary,  replacing  X  by  y  and  z  successively,  we  can  demonstrate  that 
the  same  quantity  13  independent  of  y  and  z .  Thus  it  depends  only  on  t  and  we 
can  call  it  H  (t).  If  it  is  assumed  that  H  (t)  is  continuous  in  the  domain  of 
variation  of  t,  then' there  exists  an  integral  rj  (t)  *  j  H  (t)  dt  so  that 


and  we  may  write 


_  iff  aA2  +  *  +  /  aiY 

2  1  lex/  \^dyj  [dzj 


M 

Now  (hO)  we  can  replace  everywhere  £  by  §  -  T)  (t)  since  TJ  (t)  depends  only 

on  t  and  not  on  X,  y,  t.  We  obtain 
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or,  replacing 
/ 


by  u,  v,  w^we  have 


Q+|f'-^[^+vl+w\N° 


This  equation,  together  with  the  continuity  equation  (31)  and  the  equation  of  , 
state  in  the  form  (3)  gives  three  equations  to  determine  the  three  unknowns  <|>  , 
P  and<S"in  terms  of  x,  y,  z,  t  and  arbitrary  constants. 


17 
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8.  Wave  Equation 


If  the  velocity  of  the  particles  of  the  continuous  fluid  is  small  (a  con¬ 
dition  always  satisfied  if  we  put  aside  shock  waves  of  great  intensity  or 
explosive  impulses)  we  can  neglect  the  quadratic  terms  ua  +  v*  +  w®  in  (Ul *) • 
The  latter  reduces  to 


1*2)  Q  + 


r 

-  =  0  or  V- 


dP  +  = 


a  at 


owing  to  the  definition  (35)  of  Q.  Suppose  now  that  there  are  no  external 
volume  forces;  i.e.,  V  a  0.  Then  (1*2)  becomes 


f  —  -  ai* 

J  a  at 


Now,  for  waves  of  weak  amplitude  the  mass  density  (f  is  very  nearly  a  constant 
and  for  a  variation  of  pressure  from  PQ  to  P  we  may  write 

hh)  If-  k  (p-pc) 

ro 

where  a0  is  the  average  constant  density.  Equation  (h3)  then  becomes 

a?' 

1*5)  P-P  *=  o  =_£. 

ro  uo  at 

consider  now  the  continuity  equation  (31)  where  u,  v  and  w  have  been  replaced  by 
their  expression  (3U)  in  terms  of  $ «.  Since 

.  a£'  Si .  at'  ai=  al' 

0x  3x  6y  dy  ’  3z  3z 
we  may  substitute  $  for  .  We  obtain 


1  at 


s!f  +  &  + 

dx  x  8y  ®  dz 2 


jrf 

Since  3-^1* V*  H  where  is  the  Laplacian  operator,  we  may  write  (1*6 ) 

in  the  form  V^fy^+  ^az^ ii101’*  since  0  -  oq  independent  of  x,  y  and 


Or,  since  0  -  oq  independent  of  x,  y  and  a 


W  a~° =  °«V*§' 
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According  to  the  equation  of  state  (3) ,<3^  depends  explicity  on  P.  Consequently, 
we  can  write 

^  aP 


48) 


at  ap  at 


Differentiating  (45)  with  respect  to  t  and  noting  that 6^ and  Pfl  are  constant  with 
respect  to  t  and  that  may  be  written^: 


49) 


P 


d>  \ 


O’  L- 


-  (S', 


o 


c>± ' 


Substituting  this  value  for  in  (48)  we  get 


;oi 


_  , 

at  ap 

Substituting  (50)  in' (47)  we  have  ;  . 


51) 


a 


«-0  ££•»•  fi0va§ 
”  °  at* 


a  P 


Ap 

On  setting c  =  -  equation  (51)  becomes 


52)  -j_ 

V  §  = 


i  e*i 


dil 


Equation  (52)  is  called  the  wave  equation.  It  represents  waves  propagating  in 
the  (x,  y,  z)  space  with  velocity  c,  It  can  also  be  deduced  in  the  following 
way:  suppose  that  the  volume  V  containing  a  certain  fixed  mass  of  fluid  varies 
and  becomes  V*7.  This  defines  the  cubic  dilation  £  as  the  ratio 


53) 


5- 


v'-v- 

V 


So  that  V , =  V Cl  +$ ) •  On  the  other  hand,  if  the  density  of  the  part  of  fluid 
contained  in  a  fixed  volume  V  changes  so  that  the  density^ becomes  the 

condensation  £  is  defined  by 


54) 


€  = 


cr  -  <sr 

<_T 


So  that^^(H-€),  In  general,  changes  in  density  and  volume  in  fluids 
result  from  pressure  changes.  The  excess  pressure  P-P0  over  the  static 
pressure  PQ  is  called  the  acoustic  pressure.  The  bulk  modulus  B  is  defined  as 


55) 


B=  - 


i  ( 


AsV^-V^tnbe cause  of  the  conservation  of  mass,  using  (53)  and  (54)  w®  get 

.W>1  (14-r)  (|+  eV^(u54€+S€.).Neglecting  the  cross-product $€  we 
^  ;  J  and  therefore  B  »J2,  Next, 

19  € 


-m 


I  4-S-+-  €  ^  I  oV  €--S 
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consider  a  surface  3  and  a  volume  element  Pdz  and  suppose  for  simplicity  that 
the  fluid  is  deformed  only  in  the  z  direction.  A  particle  of  the  fluid  which 
was  at  z  at  time  t  undergoes  a  displacement  dz  =  dl.  The  displacement  of  the 
particle  occupying  the  position  z  +  dz  at  some  time  t  will  be  1  +  dl  where  dl  = 
[  3z  i  Z*  *  reSfure  z  ProSPurR  z  +  dz  is  P  +  dP  - 

P  +^^jdz,  where  P  is  the  force  per  surface  unit.  Consequently,  the  net 


force  in  the  z  direction  is  equal  to  PP  -  S  (! 
mass  of  fluid  contained  in  the  volume  clement 
dynamic  equation  in  the  z-  direction  is 


(P  ♦(fgW)  -  -  S 
t  £dz  is<£Sdz.  Thus 


dz.  The 
Newtonian 


■is) 


dz  =  f^dz  3— ^ 

at* 


hence  -  =  o  .  Now,  since  P  E  B€,  we  .deduce  by  differentiation,  B  being 

a  constant:  ^  J 


iE  =  n  ^  =  -  b 

a  a  •  0z  "  ~3z 


vy  -V 


=  To  that 

O  Z 


3^  - 


B  and  we  have 

az* 


aS  =  /  e*i 
m  U  j  at* 


which  is  a  wave  equation  in  the  z  direction,  emphasizing  the  role  of  the  coefficient 
B.  For  adiabatic  processes  in  ideal  gases  we  have  Pv®*’  =  constant  so  that  we  can 
define  the  adiabatic  bulk  modulus  by  differentiation: 

OC  .  Ct*» 

v  dP  +  ct,  v  Riv  =  0. 

Thus  dP  =  -  cC  ^  P  =  -d,f  P  si  tha+.®(P  =  -  ^  =  B  and  c  =  '  a  ®LE ,  We  Can  verify 

v  J  <S 

that  (£6.)  or  (5?)  represents  a  wave  equation  in  the  z  direction.  In  that  case, 

N  the  general  solution  to  this  equation  is  h  *  f  (z  +  ct)+  g  (z-ct)  *  f  (u)  +  g(v) 
on  setting  z  +  ct  =  u,  z  -  ct  =  v.  Differentiating  h  with  respect  to  t  and 
z  we  obtain: 

§ii=5_i!3u  +  4s^z  =  cf  —  - 

at  3u  a  t  a  vat..  ^au  av  y 

£?h  =  cr^a_u  1 

at  *  1  3ua  at  dv*  at  j 
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The  functions  f  and  g  represent  any  differentiable  function.  According  to 
the  equation  h  =  f  (z  +  ct)  the  deformation  f  (z)  which  at  tirAe  t  =  0  was  at 
z  *  zo  is  at  time  t  at  z  +  ct  with  the  same  functional  form  f  (z)  ,  Thus  the 
wave  propagated  from  z  to  z  4  ct  with  velocity  c.  In  the  same  way  a  wave  with 
the  functional  fo^m  g  (z)  is  propagated,  from  z  =  z.  to  z-ct  with  velocity  c  in 
time  t.  In  the  general  case,;  consider  the  wave  equation  (52)  where  c  is  a  con¬ 
stant  and  we  take  as  a  single  independent  variable^*  =  0(  x  +P*y  +Yz  +  ct  where 

its.  The  wave  eauation  hpcnCi:  v 


.  - - -  a  L.  .UIIOJIC  u  • 

^ ^  anc^  Y  sre.  constants  .  The  wave  equation  becomes 


2. 

Thus,  when  ot  + 
(52),  Also,  Yt  = 
$  =  f ,  9U0  >f 

functions  - 


Vex1 


^1>  4- 


=  O 


=  CXx 


Y  =1,  any  differentiable  function  §  is  a 
v£y  z  -  ct  leads  to  a  solution  of  (52). 


solution  of 
Thus, 


2.  is  a  solution  where  f,  ,  are  arbitrary  differentiable 
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9,  Wave  Equation  in  Cylindrical  Coordinates 

For  our  aim,  which  is  to  analyse  the  behavior  of  sound  waves  in  tubes,  we 
have  to  express  the  wave  equation  (52)  in  a  cylindrical  coordinate  system  instead 
of  the  usual  rectangular  system.  The  cylindrical  system  is  more  adapted  to  our 
problem  and  has  the  supplementary  advantage  that  we  can  define  in  it  separated 
solutions- i, e, ,  solutions  formed  by  a  product  of  functions,  each  factor  of  which 
depends  only  on  one  coordinate.  These  separated  solutions  satisfy  a  system  of 
ordinary  differential  equations  -  which  are  much  easier  to  handle  then  partial 
differential  equations  -  and  all  solutions  of  the  partial  differential  equation 
can  be  determined  as  linear  combinations  of  each  term  of  the  set  of  separated 
solutions.  Thus,  to  integrate  equation  (52)  is  equivalent  to  determining  the 
separated  solutions.  The  coordinates  for  which  such  solutions  exist  are  called 
separable  coordinates. 

We  choose  as  coordinates  a  directed  radius  r  smarting  from  the  origin  of  the 
rectangular  coordinate  system  and  situated  in  the  (x,  y)  plane  and  angle  Q 
which  is  the  positive  angle  determined  by  r  with  the  x-axis.  The  (x,  y)  plane 
can  be  regarded  as  defining  a  cross-section  of  a  certain  cylindrical  tube.  The 
z  axis  in  this  cace  is  along  the  longitudinal  axis  of  the  tube. 

To  express  (52)  in  cylindrical  coordinates  we  notice  that  r, 0  are  simply 
polar  coordinates  in  the  (x,  y)  plane  su  that  the  definition  of  the  cylindrical 
coordinates  is  given  by 

57)  x  =  r  cos  G 

y  «  r  sin  G 

z  =  z 

or,  in  terms  of  x  and  y 

58)  V'-  \ J 

©  -  ±osr\ 

Differentiating  (58)  we  get 

% 

_  ^  —  =  Sl/A  © 

'  V  y  5 

'i  =  __L  SuaQ 

=  ~  'V-a+^a  ~  ^ 


c>lr 

& r 


=  Cos  O 


Av  / 


A/ 

/V 


•  J 


Cos  © 
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The  wave  function^  (x,  y,  z,  t)  when  expressed  by  means  of  the  variables 
(r,  ©  ,  z,  t)  becomes  another  function  ^'(r,0#  z,  t) .  In  order  not  to  needlessly 
complicate  the  notations  we  shall  retain  c|>  in  any  case.  Applying  the  formulas  of 
partial  differentiation  we  get 

6°)  _5f__  >  |  je 

6%  ~  ay  a% 

«g>Q 

©y  "  av*  ay  ay 

Substituting  (59)  into  (60)  we  get 


61) 


Cos  ©  - 


ay 

Consequently 


63) 


su©  + 

©r 

r 


ai  Sty.  9 


a© 

<3© 


cS)^2- 


a. 

djL. 

JL 


ay- 

©'r 


Cos  ©  - 


Vs 

QjQS  © 

\r 

51 


Sov\  © 


ay^ 

differentiating  (63)  we  have 
64)  c>a3> 


SiAA©^ 


a(fS  ©v-  j_  ©c 

°l  r  ^ 
+d@Lair 


Cbs© 


L 


a\rx 


Cos  ©  + 


a© 

r  1 

©§  Cos  ©  | 

a© 

V* 

ScAA© 

fir 

a© 

Y~ 

a%> 

si. 

SUa©  7 

©0 

Vs 

J 

51 

S©'© 

aai 

a© 

©v^a© 

S^e 

V" 


Cos  G 


Suvv© 

v- 


ocfo  ^ 

— — :r-  Cos  ©  ~  T —  5cv\  ©  - 
©r©0  ©r 


\2t 
c9  CD 

©e" 


Su^0  ©$  C^©’ 


a©  c 


a2! 


F  SiM.  c;  •  ft  , 

a'ya-=drj_5r  ^  +  dS 


Cos©  j^r  ,  c a  ^5  Cos©]  5© 

r: - r~  k~:  "K  J  3tx  ^  w  9 +- -^r—  -r~ 


r  Jay  a©l_d^ 


ae  r 
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.[■ 

sin 

8. 

3  ,  a  f  cos  9 

djf  cos  9  j 

L 

a-1 

093r  r 

o9 

r  1 

cos  9 

f  5-^  sin  9  +  £#cos 

e  ♦  f!f. 

cos  9 

r 

I. 0r09 

dr 

09* 

r 

Adding 

0x* 

‘  a«# 

and  —a  we 

ay52 

get 

65) 

ai*  + 

+  i  2# 

+  L 

*0X® 

0y*  dr"* 

r  dr 

r* 

sin  9  + 
dJF  sin  9 

‘  55  “ 


so  that  the  wave  equation  (5 2 )  takes  the  form 

a!#  +  i  a#  +  l.  f!£  ♦  . 

3r*  r  0r  r*  39 *  3z* 


66) 


1 


ci* 

5T> 


i 
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10.  Separated  Solution  of  the  Wave  Equation^ 

We  have  already  seen  that  in  cylindrical  coordinates  equation  (66)  is  separable, 
consequently,  a  solution  to  this  equation  is  given  by 

67)  <j)  (r,6,  z,  t)  =  A(r)  B(0)  C(z)  D(t) 

Each  factor  of  the  product  being  related  to  one  variable.  Substituting  in  (66) 
the  product  (67)  we  obtain 

a>2A  aA  ,  .  a2- b  .  ,  a2c 

68}  BCD  f  BCD  +  ACD  ^  ^  ABO 

=  ABC 

Dividing  every  term  by  ABCD  (68)  becomes 

i  s2a  ,  _u  4.1 J_  <£n 

695  A  dr2  Ar  Sv-  Bt"’-  a©2  c  az2  "D C.a  at2 

The  left-hand  side  of  (69)  does  not  depend  on  t  since  it  is  only  formed  by  the 
functions  r,0,  z  and  their  derivatives.  On  the  contrary,  the  right  hand  side  of 
(69)  is  only  formed  by  a  function  of  t  and  its  second  derivative.  Thus,  when  t 
varies,  the  left  hand  side  remains  constant  and  when  r,Q,  z  vary  the  right  hand 
side  remains  constant.  As  the  right  hand  side  is  equal  to  the  left  hand  side  we 
conclude  that  it  is  equal  to  the  same  constant,  that  we  can  term  k/*-.  Hence 


70) 


Dca  at2' 


or 


<£d 

at2- = 


D 


The  general  solution  of  this  ordinary  differential  equatior  of  the  second 
order  with  constant  coefficients  is  of  the  form 


71) 


.  kct  i  -kct 

D=  h,e  +h2e 


where h.andln^  are  arbitrary  constants.  We  see  that  (71)  represents  harmonic 
waves  only  when  k  is  purely  imaginary.  Hence  we  may  set 


72)  k&  =  eTTjr  OV-  k=  avjv/c  =  2-Trj  /'X 


where  A  is 
general,  k 
becomes 


the  wave  length  of  the  periodic  phenomenon  andV  is  the  frequency.  In 
is  complex,  k  =  a+jb,  and  involves  an  attenuation  factor  "a”  so  that  (71) 
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a  and  b  being  real  numbers.  The  arbitrary  constants  hx  and  ha  in  general 
complex  numbers.  Putting  the  right-hand  side  of  (69)  equal  to  Iv*,  we  can  write 

73)  1  2£A  +  L  3A  +  _1_  8SB  „  .  1 

,  A  0r®  Ar  Sr  Br*80fl  C  8z® 

We  find  ourselves  in  the  same  situation  as  previously.  The  right-hand  fide  of 
(73)  is  only  a  function  of  z  and  the  left-hand  side  only  a  function  of  r  and  Q. 

As  they. are  equal  to  each  other,  they  can  be  equated  to  the  same  constant (k')£ 
Thus 

7i»)  .  (k*  -(k)*)c 

The  general  solution  of  this  ordinary  differential  equation  of  the  second  order  is 
7?)  C  =  Cl[e  ' 

Replacing  the  right-hand  side  of  (73)  byfc7)eand  subtracting  - ,  we  obtain 

Br^9e 

76)  'i  £a  +  l  3a  jLa!e 

A  3rfi  Ar  3r  Br*  8©* 

multiplying  by  re 

77)  ra/l  +  L  §A  -  1  §5 

l  A  3r*  Ar  dr  K  *  j  B  80® 

Again,  the  right-hand  side  of  (77)  is  only  a  function  of  9  and  the  left-hand 
side  is  only  a  function  of  r.  Thus,  they  are  both  equal  to  the  same  constant 
that  is  customarily  called  Accordingly,  we  have 

101 

The  general  solution  of  which  is 
79)  B  =  Bi  e  +  B,,  e 

Bi  and  Be  being  arbitrary  constants,  m  must  be  a  real  number,  otherwise  (79) 
would  give  vanishing  or  unbounded  solutions.  If  m  is  not  an  integer, the  wave 
B^ejro®  +  Bae-lrc®  represents  a  progressive  wave  in  the  9-space.  For,  B  is 
changed  to  B1#  B  by  a  2K  rotation.  When  m  is  an  integer,  for  8  =  0  or  n  2n, 
n  an  integer,  b1  -  Bj.  +  Bx  =  B  remains  the  same, Physically ,  this  means  that  the 
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wave  front  after  a  complete  rotation  of  n  2 ft  degrees  comes  back  to  the  same  phase. 
The  wave  is  thus  rtationar^J’in  the  9-space, 

1 2 

^Next,  replacing  in  (77)  the  right-hand  side  by  m2,  substituting  for(k)  ,  h*  = 
Ok ’/and  conveniently  arranging  the  terms  we  get 


80) 


d*A  , 
dre 


l  ^  +  (h 

r  dr 


Finally,  replacing  r  by  s  =  hr  we  obtain  the  general  Bessel  equation 


81)  ~  +  (1-  HL)  a  *  0 

dss  a  ds  s* 

The  partial  derivative  3  /d  s  has  been  replaced  by  the  total  derivative  d/ds 
since  s  is  the  only  variable.  The  general  solution  of  (8l)  can  be  given  in  terms 
of  the  Hankel  functions  of  order  m  of  the  first  and  second  kind  ,  Hm  (2): 


82)  A  =aOlm(1)  (s)  (s). 

In  practical  applications  the  Hankel  Functions  are  used  to  represent  progressing, 
cylindrical  and  spherical  waves  for  the  following  reason.  If  we  assume  that  s  is 
very  large  so  that  we  can  neglect  l/sR,  (8l)  is  simplified.  On  setting  A sfs  =  f, 
it  becomes  in  this  approximation 

83)  li  +  f  0  0 

ds* 

This  equation  has  the  particular  solutionsOCe  representing  progressive  waves  in 
the  s-direction.  In  fact,  for  large  s,  the  Hankel  functions  have  the  asymptotic 
form 


j^C1)  =  eJ  (s  "  h  m7l_  471) 

11.  (2)  ._/T“  e-j  (s  -  ’  mjt  -  Iv.) 
1  m  V  ^  e 


-ft<arg  s<^ 

-  2ft<arg  8<n. 


Hence,  they  differ  from  f/y^T  only  by  a  constant.  Now,  we  can  define  the  Bessel 
functions  Jm  (s)  by  the  relation 


8h) 


(s)  +  H 


(2) 


m 


(s) 


27 


rw* v 


lab.  Project  FR-68 
Progress  Report  II 
Enclosure (1) 

Because  of  the  linearity  of  equation  (81)  they  are  also  particular  solutions  of 
this  equation*  According  to  (3h)  their  asymptotic  form  for  large  s  becomes 


8S)  Jir.(s)  *  i 


wtr  TT| 


.12  /  mTl  H  \ 

»Af~  cos  (s - -  -  ) 

2  hJ 

(8$)  is  more  adapted  to  the  description  of  stationary  waves  as  it  involves 
waves  in  the  +,  a  directions,  VJhen  m  is  not  an  integer,  the  general  solution  of 
(81)  takes  the  form 


86)  A  -  ^Jm  (s)  +  (s) 


where 


(*) 


(f)  2 ; 

yi*o 


/  *\?n  n 

(?)  (-1) 


r  r(w-fn+i) 


'-<*>  ■  ( » TZ  r, 

where  F  is  the  gamma  function. 


(if  (-D 


T (-m+n+1) 


When  m  is  an  integer,  J_m  (a)  degenerates  to  the  form  (-l,1"  Jm  (e),  making  Jm  (e) 
and  .J.m  (s)  linearly  dependent.  As  we  must  hnve  two  independent  solutions,  the 
second  solution  is  now  given  by  the  Neumann  function  Np,  (s),  defined  in  terms  of 
the  Hankel  function  by 


%  (®)  *  *rj  ^00  "  nTn 

■  ln  f.)  4  (»> 


1  (m-n-l)l/s\2n-m 


*  *  nl  U 

n*o 


- 1%  (.1)n  fi-USs'  (A?.  .  .  . 

11  n-o  n|  (n+m)l  \v 
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1,  2,  3  .  . 


0.^7721 


l+"  +  “+.  •  i  +  *  j  ho 

J  J  k 


so  that 

87)  A  =o(1  Jm  (s)  +  0(2  N,n  (s) 

The  Asymptotic  form  of  Nm  (s)  is  proportional  to  sin  (r-  ^  -  JL)  so  that  the 
Neumann  function  represents  a  generalization  of  the  sine  function.  1  When  h^  is 
negative  nothing  important  is  changed  in  the  form  of  the  solution,  although  the' 
Bessel  function  depends  on  the  imaginary  argument  of  s,  hr  =  ,jk1r .  For.  a  solution 
of  (81)  is  also  given  by 

.inn 

88)  Im  (.Is)  a  Jm  ( jk*r) 

where  Im  (js)  is  a  real  function  of  k*r  which,  when  m  it  an  integer,  degenerates  to 
(±  1,  ±  j)  Jm  (s). 

The  Bessel  equation  (8l)  is  solved  by  a  series  expansion.  On  setting  A  =  saC]  + 
j^a+lCp  +  ....+  s'3+nCn+i  +  .  .  .  .  ,  differentiating  A  and  equatinr  the  coeffi¬ 
cients  of  the  same  power  of  s  we  obtain  for  the  Bessel  function  of  m  order  (m  an 
integer)  the  series 


89)  jm(s)  ,  L  (|)" 
ml  <■ 


i  m+2 

— .  (£)  + 
(m+1 )i  V  * 


kl.(m+k)l  2 


s  m+2k 
(s)  +  • 


Series  (89)  converges  for  finite  values  of  s,  oscillating  as  the  trigonometric 
functions  do,  but  with  an  amplitude  which  decreases  as  s  increases.  Between  the 
first  two  Bessel  functions  J0  and  Jj  we  have  the  relation 

90)  /ja  do  =  -  J0  (a) 

d 

Remembering  that  the  fundamental  particular  solution  of  (66)  is  (67),  we  have  to 
substitute  for  B,  C,  D  the  values  (71),  (75)>  (79).  For  the  radial  function, 
if  we  consider  the  progressing  waves  we  can  choose  (82). 

As  we  need  only  two  constants  we  write  tljo  particular  solution  in  the  form 


(91)  %  ■  (^Hm(1 )  *  0^))  e±jb0t  * 
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If  m  is  not  an  integer  this  describes  a  progressing  wave  that  propagates 
in-  the  z  direction  in  a  screw-like  manner.  VJhen  m  is  ah  Integer  we  may  also 
represent#  by  a  superposition  of  stationary  waves.  In  this  case  we  choose 


for  the  radial  wave  function  the  solution  (87) ,  so  that, 

92)  f,  -  (<V„  *  **  Hm)  e^me  e^bot  #')*-**"  z 


The  general  solution  is  a  linear  combination  of  (91)  or  (92)  and  may  be 
written  for  instance  using  (92) 

*  . 

93)  i-  y  (V»  ♦«.%)  •***  a4*"*  e*  ^  *  2 

&tp  k* 

1 

m,  k  ,  k  being  variable  parameters. 
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11,  Digression  about  Function  fpaces 


6‘ 


It  would  appear  surprising  that  the  solutions  of  certain  differential 
equations  could  bo  expressed  in  several  different  ways,  for  instance  that  the 
general  wave  solution  of  (8l)  could  be  expressed  either  in  terms  of  elementary 
progressive  waves  or  in  terms  of  stationary  waves,  Duch  a  property  follows 
from  the  fact  that  the  elementary  functions,  for  instance  Jm  or  Nm  or  in  the 
case  of  a  Fourier  series  a0,  an  cos  nx,  b^  sin  nx,  form  what  is  called  a  com- 
w  plete  system  of  functions.  This  is  readily  understoou  by  apjealing  to  geometrical 
intuition.  Consider  a  linear  vector  space  of  n  dimensions,  Rn.  This  means  that 

.  .  . ,  en  such  that 


n,  for  instance  p  can  be  represented  in  terms  of  the  or 


there  exists  a  fundamental  system  of  unit  vectors  e  ,  eP, 
all  other  vectors  of  T  * 

more  precisely  can  be  regarded  as  linear  functions  of  e^: 

_  _  .  %  _ 

9h)  P°pe+p2ea+...+  pnen  =2,  ^k  ek 

■fe  » \ 

We  say  tha-  the  e^  form  a  complete^ system  of  basic  vectors.  Of  course,  we 
can  choose  another  set  of  unit  vectors  e^  which  will  also  be  complete.  'This 
means  that  the  "p^  can  be  expressed  as  linear  functions  of  e^; 


9$) 


)  t 


7  7 


p  =  pi«i  +  p^Gp  + 


/  / 
+  pnen 


'VO 


/  / 
^kek 


asso'’ 

—  t 


We  can  also  say  that  the  set  p  has  the  components  p'^  in  the  representation 
-iated  wittL  t.he*e"'k  and  pk  in  the  representation  asfociated  with  e^  or  that 
e  k  define  two  different  coordinate  systems.  The"p,  and'eyK  are  quantities 
that  are  defined  by  a  magnitude  and  a  direction.  There  are,  however,  mathematical 
objects  of  a  more  general  character  which  satisfy  the  same  conditions  of  completeness 
and  linearity  and  consequently  define  a  linear  space.  The  elements  of  this  linear 
space  instead  of  being  vectors  are  functions.  For  instance,  we  can  speak  of  tne 
space  of  continuous  functions  f(x)  for  0<x<l.  If' we  take  as  basic  vectors  the 


quantities  1,  cos  x,  cos  2x, 


cos  nx 


sin  x,  sin  2x, 


sin  nx,  .  .  ,  ,  this  system  of  functions  is  complete  with  respect  to  uniform  con¬ 
vergence  in  the  space  of  continuous  periodic  functions.  This  means  that  any  con¬ 
tinuous  periodic  function  f(x)  being  given  in  a  certain  interval,  it  will  always 
be  possible  to  represent  such  a  function  in  the  following  way: 


09 


•O 


96) 


f(x)  =  ao  +  2  ak  cos  104  +  ^  ^k  s^n  k* 

■feWi  fe-i 

We  recognize  the  expansion  (96)  as  a  fourier  series.  But,  on  the  other  hand, 
the  system  of  functions  1,  x,  x2,  .  .  .  ,  x11,  ...  is  complete  with  respect  to 
ordinary  point-wise  convergence  in  the  space’of  continuous  functions.  Thus,  if 
in  example  (96)  f(x)  represents  a  continuous  periodic  function,  it  can  also  be 
represented  in  terms  of  the  basic  functions  1,  x,  x^,  .  .  .  . 

These  examples  are  sufficient  to  show  that  under  very  general  conditions  a 
continuous  function  can  be  linearly  represented  in  terms  of  a  sequence  of 
definite  other  functions  and  that  the  choice  of  these  functions  could  not -be  unique. 


31 


k  Lab.  Project  FR-68 
‘  Progress  Report  #1 
Enclosure  (1) 

For  instance,  a  function  could  be  represented  in  terms  of  the  Legendre  polynomials, 
or  in  terms  of  the  Hermite  polynomials.  More  particularly,  the  Bessel  function 
Jm(s)  where  tn  is  an  integer  form  a  complete  system  able  to  represent  under  general 
conditions  a  continous  function  f(s).  A  general  wave  solution  of  (8l)  can  thus 
be  represented  in  terms  of  progressive  waves  or  in  terms  of  stationary  waves. 
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7 

-  1? .  Proper  Frequencies  of  Radial  V/aves 

Mo  simple  result  could  be  obtained  from  the  study  of  waves  as  complicated 
as  those  described  by  (93)  that  have  variable  frequencies  and  amplitudes. 

Experiments  to  test  the  acoustic  properties  of  a  sample  at  the  termination  of 
a  tube  are  necessary  using  monochromatic  plane  waves,  like 

97)  f (a,  t)  =  Ae3kz  e3kct 

with  constant  amplitude  and  frequency.  In  this  way,  we  are  able  to  unequiv¬ 
ocally  compare  the  amplitude  and  phase  after  being  reflected  by  the  sample  with 
the  incident  ones.  Thus  the  experimental  problem  is  to  produce  only  waves  of 
the  form  (97)  and  not  of  the  general  form  (93).  This  problem  can  be  discussed 
mathematically  because  (97)  is  a  degenerate  form  of  (93).  Thup,  we  have  to 
determine  under  which  conditions  (93)  reduces  to  (97).  Thysically  production  of 
waves  like  (93)  means  that  it  is  not  sufficient  for  producing  a  monochromatic 
wave  to  give  a  piston  at  the  end  of  a  tube  a  sinusoidal  motion  u  =  Uo  sin  2nv*t 
because  Huyghen's  principal  states  that  each  point  of  the  piston  surface  can  be 
regarded  as  the  source  of  a  wavelet  expanding  in  the  three  dimensions  of  space 
and  thus  giving  rise  in  general  not  only  to  waves  along  the  z  axis  (which  is 
the  tube  axis)  but  to  waves  in  the  perpendicular  (x,  y)  plane  that  are  described 
by  Bessel's  function  and  consequently  do  not  have  constant  amplitude  and 
frequency.  This  part  of  the  general  wave  is  just  what  is  described  by 
°Cim  Jm  +c^-em  Nm  i-n  (93).  As  s  =  hr  it  is  in  the  direction  of  the  polar  axis 

r  «=Vx*  +  y*  which  is  perpendicular  to  the  z  axis. 

Now  these  radial  waves  are  associated  with  a  discontinuous  series  of 
freouencies,  because  they  have  to  satisfy  the  two  following  conditions  that  are 
very  restrictive: 

1.  Their  amplitude  cannot  be  infinite  at  the  origin  r  =  0,  which  would  be 
contrary  to  experimental  results. 

2.  They  must  be  reflected  by  the  walls  of  the  tube,  regarded  as  a  perfect 
reflector. 

The  first  condition  eliminates  the  Neumann  functions  Nm(s)  which  are  not  finite 
at  r  s  0. 

The  second  condition  requires  that  the  components  of  velocity  of  the  radial 
waves  at  the  wall  of  the  tube  vanish  so  that  the  waves  are  entirely  reflected. 

Thus,  if  the  radius  of  the  tube  is  defined  by  r  =  a,  we  have,  as  the  velocity 
components  derive  from  a  potential  $  according  to  (3h) 


■  "(If)*  ■  T*  =  -  (§)a  =  °* 
Hence,  by  solving  (6l)  we  obtain 
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According  to  the  first  condition  we  have  to  put  in  equation  (93)  Nm  -  0  for 
all  m.  Further,  taking  the  derivative  of  (93)  with  respect  to  r  and  letting 

=  0,  we  obtain 

98)  =  0  for  r  =  a  or  for  s  =  ha 

\  3r/r  =  a 

a  being  givan,  corrixiion  (98)  cannot  be  satisfied  for  any  arbitrary  value  of  h  or  s, 
but  only  for  a  sequence  of  discontinuous  values  hn  or  sn.  Thus,  if  we  term  hm  the 
parameter  associated  with  the  function  of  m  order  Jm,  the  different 'values  of  hm  for 
which  (98)  is  satisfied  can  be  denoted  hmn,  corresponding  to  smn.  In  the  following 
table  we  give  the  values  of  8mn  corresponding  to  the  zeros  of 

4  for  m  =  0,  1,  2,  3 

and  n  =  0,  1,  2,  3# 


TABLE  1 

n 

<M 

»xn  for  Jo=0 

s2n  for  =0 

s,n  for  j;  .0 

3*n  for  Jl=0 

0 

0 

l.Slj 

3.05 

lj  .20 

1 

3.83 

8'. 33 

6.70 

8.02 

2 

7.02 

8.5h 

9.97 

11.3b 

3 

10.17 

11.71 

13.17 

.  lU. 59 

Consider  the  general  solution 


Qr\  ^  cy  ,  ,h  ^  ^m9  ±jbct  z 

93  )  ¥  Z,  Jm  (hmnr)  °  e  e 

m  k,k‘ 

In  order  that-  $  be  independent  of  r,  which  is  the  first  condition  for  a 
plane  wave  along  the  z  axis,  hmr)r  must  reduce  to  a  constant,  otherwise  the  function 
Jm  (hmnr)  cann0't  reduce  to  a  constant.  The  term  hmnr=  C  is  impossible  since  hmn 
and  r  vary  independently  of  each  other,  unless  C  =  0,  which  implies  that  hmn  -  0 
since  r  always  varies  from  r  =  0  to  r  =  a.  Thus  =  0.  Referring  to  the  defini¬ 
tion  of  J-  given  by  (89)  we  see  that  this  solution  is  satisfied  only  by  Jm  =  0  when 
m^O,  so  that  in  this  case  wave  ( 93 ^ )  vanishes.  For  m  =  0  we  have  the  solution 
Jc  =  1  for  any  value  of  r  and  also 

^0  =  r  !  =  0 

'  3r  ° 

Thus  all  conditions  are  satisfied.  As  h  2  =  -(k k7  **  0,  (93)  reduces  to 


Jkz 


99)  £ 

i.e,,  to  the  monochromatic  plane  wave  (97).  This  lasts  until  we  meet  with  the 
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least  value  of'  which  is  s10  corresponding  to  J-,  =  0  according  to  Table  1. 

To  explain  this,  consider  pure  radial  waves  along' the  r  axis  so  that  there  is 
no  preparation  along  the  z  axis.  According  to  (75)  this  means  that  k  -(k  ;  - 
Taking  account  of  (7?)  and  (71)  we  see  that  h  corresponds  to  the  frequency  of  a 
periodic  phenomenon  f  such  that 


100) 


V*  = 


he 

271 


so  that  D  =  h-^  +,hp& 


how  Table  1  shows  us  that  the  smallest  value  of  8^  =  h^a  is  at  Spo  =  1.8U» 
Thus,  the  smallest  value  of'  h^  is  h^Q  =  l.Sb/a  and  the  smallest  radial  frequency 

is 


101) 


t  _  1.8lic 

^0 =  1ST 


Similarly  the  second  possible  frequency  V* corresponds* to  s2o  in  Table  1  so  that 

-  3.05c 

^20  2% 

3  83c 

is  the  following  frequency.  Similarly,  the  next  admissable  frequency  is  Vbl  =  ^ - 

Tn  general,  the  different  radial  frequencies  are  obtained  by 

v-mn  -  hmr,o/2’i  =  »mlo/2«a. 


Thus,  we  have  to  stay  away  from  these  frequencl  ?s  if  we  want  to  operate  on  a  pura 
monochromatic  olane  wave; 
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13 •  Definition  of  Specific  Acoustic  Impedance  -  Comparison  Between  Hater  and 

Air  Filled  Tubes?* 


The  wave  equation  (£2)  is.  satisfied  by  potential  <§  and  also,  as  equation 
(£6)  has  shown,  by  displacement  h.  Morever,  the  same  equation  is  satisfied  by 
pressure  P,  local  velocity  (u,  v,  w)  and  density  o  .  To  show  that  it  is  approxi¬ 
mately  .satisfied  by  pressure  P,  we  only  nave  to  use  equation  (b$),  where  o0  if 
a  constant  and  §>*  is  simply  called  For,  differentiating  (£2)  with  respect  to 
time  t, . multiplying  by  the  constant  o0  and  taking  account  of  the  commutativity 
of  the  differential  operators 


3_  a. 
ax-  ay’ 


—  and  4:  since  § 

az  at 


is  a  uniform  function: 


102) 


*  -  1 

V  p  &  at* 


Similarly  the  velocity  depending  on  the  potential  §  we  get  for  u 


103) 


1 

c*at */ 


or 


<L*)or 
dx  /  . 


and  similarly  for  v  and  w: 


w 


i  ii 

cz  at® 


a 

V  v 


_i_  afv 

C*  0t® 


Consider  the  simple  cafe  where  a  monochromatic  wave  proparatps  in  a  certain 
indefinite  medium  in  the  positive  z  direction  so  that 
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1*)  *(*,  t)  -  Ae^  eJkct 

from  (b5),  the  acoustic  pressure  is  given  by 

105)  P *P os  Cq  *  Oq  jAkb  e^k'2  eJkct 

and  the  velocity  w  of  elements  of  the  medium  in  the  z  direction  is  given  by 

106)  w  s  -  =  -  jAk  e^k^  e^kct  or 

3z 


jkcfcc  3z 


The  quantity  ejkct  appears  everywhere.  •  To  avoid  using  variable  quantities  of 
this  kind,  it  is  convenient  to  define  the  new  concept  of  specific  acoustic  imped¬ 
ance  of  the  medium  transmitting  the  plane  wave  as  the  ratio  of  the  external  pressure 
■P  -Po  to  the  local  velocity  w.  It  represents  the  reaction  of  the  medium  to  the 
external  pressure  and  is  a  certain  measure  characteristic  of  this  medium.  For 
instance,  if  no  velocity  is  communicated  by  external  pressure  to  the  particles 
of  the  medium  the  specific  acoustic  impedance  is  infinite.  It  is  defined  by 

107)  Z  »  «  o  c 

w 

according  to  (105)  and  (106)  for  a  plane  wave. 

He  see  that  finally,  despite  the  intervention  of  the  macroscopic  quantity  w, 
it  can  be  measured  only  by  determining  macroscopic  quantities  like  o0  and  c. 

As  the  subject  of  this  report  is  concerned  with  the  application  of  theory  to 
water  filled  tubes,  it  is  suitable  to  compare  the  accustic  properties  of  air  and 
water  as  far  as  their  behavior  in  tubes  is  concerned. 


a.  In  air  at  rest  the  sound  velocity  is  c  =  1130  ft/sec.  In  water,  c  = 
5000  ft/sec.  Thus,  the  wavelength  X  =  ^for  a  given  frequency  is  h .1*2  times 
larger  in  water  than  it  is  ift  air. 


b*  The  specific  impedance  for  air  a l  209C  is  onc  =  LI, 26b  — ^ 

r*  _  s  v  swift 

w  c 

0Qc  =  3.99xl06, 


and  for  sea  water  at  20  °C  isoQc  =  1.5x10-*  — — .  For  iron  anc^soff0  steel, 

v  j  A  m**  A  ’ 

fn  A  VIA 


cin^'-ec 


cm  -  sec 

more  than  10"’  times  the  value  for  air,  but  only  26.6  times 

Z  water, 

Z  air 


larger  than  the  specific  impedance  of  water.  On  the  other  hand,  °  wauer/  B  36I1O. 


Thus,  for  a  given  sound  velocity,  the  amplitude  of  a  pressure  wave  in  water  will  be 
36bO  times  greater  than  that  in  air. 


c.  If  the  acoustic  impedance  Zi  of  medium  1  with  local  velocity  is  much 

larger  than  the  acoustic  impedance  Z2  of  medium  2  with  local  velocity  wo  for. equal 
pressures  in  mediums  l_and  2  we  have 

1 


a.  a* 

Zi  =  7^  “  Since  zl»  z2  we  have  w2  >>wl*  Thus,  acoustic 


*2  '  Z2 
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waves  from  medium  2  will  not  penetrate  medium  1  due  to  the  smallness  of  w-, . 

Since  the  acoustic  impedance  of  steel  is  not  much  greater  than  that  of  water,  we 
can  expect  that  acoustic  waves  traveling  in  steel  tubes  filled  with  water  will 
penetrate,  the  walls  of  the  tube  to  some  degree.  The  standing  waves  actually 
induce  vibrations  in  the  tube  which  in  turn  transmits  these  vibrations  to  the 
water,  disturbing  the  plane  wave  propagation.  This  necessitates  a  study  of  the 
characteristic  frequencies  of  the  tube.  Difficulties  encountered  by  tube 
vibrations  can  be  partially  eliminated  by  increasing  the  tube  thickness*,  by 
damping  the  vibrations  by  surrounding  the  tube  with  sand  or  by  avoiding  the 
range  of  characteristic  frequencies.  The  Germans  performed  an  experimental 
study  of  their  tubes  before  using  them, 

d.  The  attenuation  of  waves  is  much  less  in  water  than  in  air,  which 
simplifies  the  equations. 

e«  The  phenomenon  of  refraction  is  important  in  water. 

f.  Noises  and  reverberation  are  obstacles  to  reception  as  in  air.  They 
are  particularly  important  in  pulse  transmission  and  reception. 

Tn  general,  the  length,  inside  diameter  and  tube  thickness  cannot  be 
arbitrarily  chosen  but  are  connected  with  the  possibility  of  definite  oscilla¬ 
tions.  The  length  of  the  tube,  in  the  case  of  standing  waves,  is  related  to 
the  maximum  wave-length  by  the  formula 

108)  L*  j 

for  a  tube  closed  at  one  end  and  open  at  the  other  or  by  L=  X/2  for  a  tube  closed 
at  both  ends.  A  minimum  frequency  of  Vmin  =  100  cps  will  correspond  to  a  maximum 
wavelength  X  s  1$2£  cm,  assuming  that  the  tube  is  infinitely  rigid  and  that  the 
velocity  of  sound  c  =  l£2,!>00  cm/sec.  Thus,  according  to  (108)  the  length  of  the 
tube  will  be  equal  to  3.8  meters. 

t 

For  a  pulse  tube  with  a  wave  train  formed  by  10  wavelengths  and  a  minimum 
frequency  =  1000  cps,  we  have  X  =  1$2,$  cm  for  each  wave  and  X  =  1^25  cm 
for  the  entire  wave  train.  The  length  of  the  tube  must  be  at  least  equal  to 
a/2  in  order  for  the  pulse  termination  not  to  interfere  with  its  front  at  the 
initial  section  of  the  tube  after  a  complete  cycle.  Other  elements  intervene 
against  lengthening  the  tube  too  much.  As  the  acoustic  impedance  of  water  is 
Close  to  that  of  steel,  prope?-  oscillation's  are  excited  in  the  tube  which  cannot 
be  avoided  by  clamping  the  tube.  Their  study  is  related  to  that  of  elasticity 
In  solids. 


♦Professor  Skudrzyk  of  Pennsylvania  State  University  has  noted  that,  as  the  walls 
of  the  tube  become  thicker,  more  of  the  acoustic  energy  of  the  sound  wave  travels 
in  the  wall,  thus  complicating  the  mechanism  of  plane  wave  propagation. 
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l)i.  Reduction  in  Velocity ^ 


The  sound  velocity  in  water  filled  tubes  is  not  quite  the  same  as  that  in 
free  water,  due  to  the  flexibility  of  the  tube.  The  sound  velocity  in  a  tube 
with  ideally  rigid  walls  being  denoted  as  c0  (same  as  in  free  water)  and  c'  oeing 
the  actual  velocity  in  the  tube,  the  ratio  c'/b 0  is  given  by  Korteweg’s  formula 


109) 


where  a  is  the  internal  radius  of  the  tube,  h  is  the  thickness  of  its. wall,  K  is 
the  Bulk  Modulus  of  Flasticity  of  water,  E  is  Young's  modulus  for  the  material 
of  the  tube.  An  improved  form  of  th^  equation  has  been  given  by  Hutte: 


110) 


/ 


c_ 

C , 


1+2 


K[(r?/n  *  3 

E[(rj/r!)?  -  fj 


1 

'2 


where  is  the  internal  radius  and  rj  the  external  radius  of  the  tube.  Even  in 
the  care  of  infinitely  large  wall  thickness,  c  /cQ  is  less  than  1,  due  to  the 
compressibility  of  the  wall  material.  For  steel,  there  is  a  limiting  value  of 
c/c0  =  0.9892.  We  can  see  that  reduction  in  velocity  can  be  corrected  by  an 
increase  of  the  tube  thickness.  For  instance,  if  rg  =  2r^,  formula  (110)  yields 


HI)  ~ 
co 

If  r2  “  \  rX 


=  0.982 

—  =  0.972. 
Co 


Thus,  the  thicker*  the  tube,  the  better  the  experiment,  as  the  reduction  in 
velocity  becomes  negligible.  One  must  also  consider  the  fact  that  some  acoustic 
energy  travels  in  the  walls  of  the  tube  and  hence  tube  thickness  will  influence 
proper  radiations  in  the  tube. 

As  far  as  the  internal  diameter  is  concerned,  .it  should  be  less 
than  the  minimum  wavelength  for  two  reasons: 

First,  the  tube  through  which  the  wave-train  proparates  can  be  regarded  as  an 
obstacle  to  this  train.  When  the  wavelength  is  very  small  compared  to  the  linear 
dimensions  of  the  obstacle  there  is  a  shadow.  When  the  linear  dimensions  of  the 
obstacle  are  of  the  order  of  the  wavelength  of  the  vibration,  diffraction  occurs, 
so  that  propagation  of  a  train  of  plane  waves  along  rays  is  not  possible .  When 


#See '  footnote , '  page  Jrff 

3 $ 
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the  wavelength  is  several  times  larger  than  the  oK'-kao1 r  /'in  this  ease  the 
internal  radius  of  the  tube),  the  obstacle  coes  not.  art  rcc  i ';b\y  perturb  the 
wave  motion. 


fecond,  we  have  already  reen  that  the  first  fremjerey  correrprnd  i  r,;  to  a 
radial  vibration  is  given  by  (101)  so  that  only  a  plane  wave  is  propagated 
below  the  lower  frequency  limit  V^q  given  by  (101).  /or  standing  wave  tubes, 
there  arc  strong  characteristic  vibrations  of  the  tube  wall.'  excited  by  the 
radial  oscillations  and  measurements  a  e  very  difficult  above  the  limit 
frequency  VpQ.  The  characteristic  vibrations  of  the  tube  include  longitudinal 
vibrations,  bending  oscillations  symmetrical  to  the  a  vis  and  two  and  three- 
dimensional  bending  oscillations.  The  resonance  frequencies  of  the  longitudinal 
vibrations  are  again  associated  with  Vessel's  functions.  mhey  are  given  by  the 
formula 

112 )  ^Lm  =  TiT. 

whore  L  is  the  length  of  the  tube.  For  waves  of  the  fir:  t  kind  m  >\  for  tne 
number  of  nodes  on  the  total  length  L.  ’or  waves  of  tne  seconc  kind,  m2  2, 

The  third  kind  of  waves  have  ?n  nodes  on  the  circumference  of  any  circular  ring 
(n>?)  and  r,-l  nodal  circles  in  the  entire  length,  the  characteristic 
frequencies  being  higher  than  those  of  the  second  kind  of  vibrations.  'rriir 
limits  the  length  that  can  be  given  to  the  tube  as,  according  to  (112),  the 
longer  the  tube,  the  lower  the  first  characteristic  froouency.  According  to 
equation  (112),  it  is  impossible  to  avoid  the  character  is.  tic  frequencies  of 
the  tube.  For  L  =  1°°  cm,  V? q  =  762.5  rps }for  L  =  2r,0  cm,  V[  q  *  3°1.3  cps, 
for  L  =  3f*n  cm,  Yjq  =  25)’. 2  rps.  To  avoid  suen  a  niaso/ice,  it  is  rocor- ended 
that  the  thickness  of  the  tube  (r?-r-|l  be  sn~b  that  tne  ratio  h  =  (r?-r-|  )/r-,  be 
in  the  range  0.2  -h^3.  ''oreover,  the  amplitudes  of  the  tube  vibration  should 
be  reduced  by  rubber  rings  and  a  filter  formed  by  two  rings,  fastened  together 
by  a  brazed  bronze  sheet, 

Next,  consider  formula  (101)  which  gives  the  lowest  1css.p1  frequency  above 
which  no  pure  {'lane  wave  propagates  through  the  water  of  the  tube.  From  (101) 
we  have 

Pa  =  X  =  0.566  X  or  X  =  3.h?n, 


where  a  is  the  radius  of  the  tube.  For  a  -  2.5  cm  and  c  =  l5?,500  cm/sec,  we 
haveXn,in  =  3Ji2  <  2.5  =  6,55  cm  and  vmx  =  S  =  17,650  cps.  Above  Yma; 

=  17,650  cps  we  have  radial  waves  superimposed  upon  the  transverse  waves  and 
the  wave  me y  no  longer  be  considered  a  plane  wave. 
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15.  Propagation  of  Plane  Waves  in  Water  Filled  Tubes^>^ 


We  begin  by  a  little  digression  about  complex  numbers.  Any  complex  number 
may  be  written  in  the  form  w  =  u  +  jv  where  u  and  v  are  real  numbers  and  j  =  /-I . 
The  numbers  u  and  v  may  be  regarded  as  perpendicular  coordinates  in  the  w  plane. 
Thus  a  complex  number  is  utilized  when  we  are  dealing  with  a  quantity  that  depends 
on  two  other  independent  quantities.  Passing  to  polar  coordinates  we  may  write 
w  =  r  cos  9,  v  =  r  sin  Q  so  that  w  takes  the  form  w  =  rei®  where  r  plays  the  rcle 
of  an  amplitude  and  9  that  of  a  phase  angle,  which  is  very  convenient  for  the 
study  of  periodic  motions.  As  sample  materials  in  water  filled  tubes  modify  both 
amplitude  and  phase  of  plane  waves,  their  action  is  most  conveniently  represented 
by  complex  variables.  By  definition,  |w|  =  r  =  /u^  +  Vz ,  and  9  =  tan  -IV, 

u 

Now,  returning  to  the  case  of  a  water  filled  tube  with  uniform  cross-sectional 
area  S,  starting  at  the  origin  of  coordinates  z  =  0,  and  extending  to  the  right 
along  the  z  axis.  The  harmonic  wave  produced  at  the  origin  by  an  audio  oscillator 
gives  rise  to  a  series  of  plane  waves  like  (104)  and  (106)  whose  characteristic 
impedance  is  given  by  (107) . 

When  the  tube  is  of  finite  length  L  and  closed  at  z  =  L,  the  plane  wave 
generated  at  z  =  0  will  be  reflected  at  the  termination  z  =  L  of  the  tube  and 
there  will  be  a  mixture  of  incident  and  reflected  waves.  The  total  pressure  at 
any  point  z  at  any  time  t  will  be  given  by 

113)  P  *  (P+  ejez/c  +  P_  e+;iez/c)  e*jet 

where  e  =  2^  is  the  circular  frequency.  P+  is  the  amplitude  of  the  incident  wave 
and  P_  the  amplitude  of  the  reflected  wave,  P_includes  any  change  of  phase  at 
the  termination  and  consequently  the  amplitudes  are  complex.  Using  the  definition 
of  pressure 


P 


at 


given  by  (45)  (where  P  is  written  *  r  Lead  of  P-Po  and  J  instead  of$'*)  and  for 
dg  the  plane  wave  equation 

114)  §>  =  (J+  §jez/c  +$>_  e+iez/c)  e--*et 


we  can  write  P  in  the  form 

115)  r  =  i  (<5oje)  eU«/=)  e1^ 

so  that 


116)  do  3%  = 


*Recall  that 
may  let  iL. 


<Ve§- ■ 

W  (t)  £  O 


P+  and 
P_ 

^<5  al? 

-v — ~  - — ■>  —  and  that  for  waves  of  weak  amplitude  we 

<)z  az 

so  that  ^  may  replace  <|  everywhere  in  our  equations. 
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In  the  same  way,  because  of  (34),  the  velocity  reduces  to  w  =  -  If.  Thus , 

d  z 

according  to  (114)  and  (106)  we  may  write 


w  =  (je/c)  (<j>+  ej£Z/c  -  $  _e^£i^c)  e  ±^' 


so  that 


118)  w  -  1  b  P 

j  e  d  z 

The  specific  acoustic  impedance  may  now  be  given  by 

119)  Z  =  P  =  c§+  e“jez/c  +$-  e+;iGz/c 

—  CJ  .  A— i  — — .  . . . i— 

§+  e"j eZ/c  e+jez/c 

From  (116)  we  have 

120)  Z  =  ±<£c  P+  e~j£2/c  +  P_  e+^£Z^c 

P+  e-jcz/c  _  p_  e+jez/c 

—  -  . 

At  z  =  L,  the  impedance  has  the  value 

121)  Z  =  *^oc  Po.  e~jgL/c  t  P.  e 

p+  e-jeL/c  .  e+jeVc 

✓ 

Obviously,  we  could  have  located  the  termination  at  z  =  0  and  the  initial  section 
at  z  *  ■  L  which  would  give  for  (120)  / 

**  *■ 

122)  Z  =  c  P+  +  P. 

°  P+  -  P_ 

—  - 

We  can  put  formulas  (121)  or  (122)  in  a  more  convenient  form  to  determine 
Z^  from  the  experimental  data.  We  define  the  pressure  ratio 

123)  P__  =  e~2W»  so  that 
P+ 

i>  =  low  P.  -  log  P-  Dividing  the  numerator  and  denominator  of  (121)  by  P+ 

- St - — - -  ,  °  \  J  S  T 

and  multiplying  them  by  e^we  get 
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12li) 


'gji-jeL/c  +  e-  (p-jeL/c) 
e(i+j£L/C  -  e" 


The  numerator  of  (12h)  is  simply  2  cosh 
2  sinh  (u-jeL/c)  so  that  (12h)  reduces  to 


((i-juL/c)  and  the  denominator  is 


125)  Z  =±o0c  coth  (n-jeL/c) 

The  parameters  o0,  £=  2T-V,  L  and  c  are  known.  p  is  determined  experi¬ 
mentally.  There  are  two  methods.  First,  the  method  of  the  pulse  tube  and 
second,  the  method  of  the  standing  wave  tube.  We  shall  first  discuss  the 
pulse  tube  method.  * 
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9  10  12 

l6.  Measurements  by  the  Pulse  Tube  Method  *  ’ 

i 

Let  p,  =  px  +  and  e'2|Il  =  r,  2 =  9  so  that  (123)  becomes 


126)  —  =  e“2^l  e  =  re-^ 

F* 

Thus,  r  represents  the  ratio  of  the  real  amplitudes  of  the  incident  and 
reflected  waves  at  z  =  L  where  the  material  sample  has  been  set.  It  is  an 
experimental  number  that  can  be  measured  by  the  methods  indicated  in  the  usual 
works.  In  the  same  way,  9  represents  the  real  change  of  phase  at  the  interface 
and  can  be  directly  measured.  Consequently  p  and  Z ^  are  known. 

In  (12ii)  we  can  divide  numerator  and  denominator  by  eP“j£L/c  an(^  denote 
e~2(p~jeL/c)  as ^ t  we  obtain 


127) 

Z  = 

l+n 
±a0c  — -i 

1-ti 

we  have 

128) 

ii 

_Z_  _  lla 
oQc  1- 3 

^—3  or,  setting  =  V' 


aoc 


Equation  (128)  represents  a  well-known  conformal  transformation  of  the 
theory  of  functions  of  a  complex  variable.  On  setting  +jv,  t)  =  x+jy,  the 

lines  u  =  u0  =  const.,  v  =  v0  =  constant  generate  a  rectangular  coordinate  system 
in  the  u,  v  plane.  In  the  x,  y  plane  these  lines  are  represented  by  two  families  of 
tangential  circles.  For,  replacing  in  (128)  r)  by  x+jy,  multiplying  and  dividing 
the  right-hand  side  of  (128)  by  (l-x)+jy  we  obtain 


129) 


so  that 


V-i 


u+jv 


130) 


r- 


u+jv  = 


(l+jy)2-*2 

(l~x)2  +y2 


,  2  2 
i-y  -x 


*5 


(l-x)2+y2  (1-x)2  +y  t- 

On  setting  u  =  constant  *  u0,  v  =  constant  =  v0,  we  obtain  from  (130) 


131) 


(x~  Sv1 


+r 


(l+«o)2 
2 


(x“1)2  +  (y"  ^  =  v? 


(a) 


(b) 
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Replacing  u0  in  (131a)  with  numerical  values  we  obtain 

TABLE  II 
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TABLE  III 


* 

Equation  of  Circle 

Center  a> 

Radius 

-ii 

(x  -  l)2  +  (y  +  =  l/l6 

(1,  -  %) 

1 

4 

-3 

(x  -  l)2  4  (y  +  1/3 )2  =  l/9 

(1,  -  1/3) 

1/3 

-2 

(x  -  l)2  +  (y  +  tO2  «  -4- 

a,  -  i) 

1 

2 

-1 

(x  -  l)2  4  (y  +  l)2  -  1 

(1,  -  1) 

1 

0 

(x  -  1)^  4  (y  -eo)2  =  00 

(l,oo) 

00 

1 

(X  -  l)2  4  (y  -  I)2  =  1 

a,  1) 

1 

2 

(x  -  l)2  +  (y  -  s  4 

a, 

1 

2 

3 

(x  -  l)2  4  (y  -  1/3)2  =  1/9 

(1,  1/3) 

1/3 

h 

(X  -  l)2  4  (y  -  i)2  =  1/16 

(1,  i) 

A 

4 

Thus,  lines  of  constant  v  =  v0=£o  in  the  U,  v  plane  transform  into  circles  with 
centers 


(1,  i— )  and  radii 
vo 


1 

vo 


tangent  t-o  the  line  y  =  0  in  the  x,  y  plane.  The  line  v0  =  0  in  the  u,  v  plane 
is  undefined  in  the  x,  y  plane.  These  relationships  are  shown  in  Figure  1. 
Charts  giving  values  of  u,  v  corresponding  to  values  of  x,  y  are  available 
commercially  and  are  known  as  fmith  charts. 

Returning  to  (125)  and  substituting  cosh  ix  =  cos  x,  sinh  ix  =  i  sin  x  we 
r.-ay  write 


^  AA  \ 

AXJ 


00c 


4A 

<5 


tan  h  }i.  (I  tar:2 


GL 


4. 

UCU1 


CL  /  O  \ 

H  (1  -  tan  h2  p.) 


tan  h 2  p.  4  ta«>2  ~ 


which  determines  u  and  v  in  terms  of  p  and 
can  also  use  the  series 


When  jjr 


and 

c2 


are  < 


we 


133) 


tan  h  (i  »  (i  -  +  2  -  17  * 

3  1<  315 


tan  & 


£  -  Jr  +  2  + 


3c 


T 


17sl7  , 

15c*  31^0? 
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£L 


Equation  (132)  in  a  first  approximation  for  high  values  of  the  impedance  with 


p,  less  than  0.1  reduces  to 
c 


13U) 


r- 


ud  *  j  fi-i1®) 


C  * 


These  formulas  give  the  characteristic  impedance  of  the  water  at  the  inter¬ 
face  z  =  L  separating  the  material  under  investigation  and  water.  It  is  also  of 
great  interest  to  know  the  value  of  the  acoustic  impedance  inside  a  layer  formed 
by  this  material  which  is  supposedly  homogeneous. 

This  has  been  done  as  follows.  Consider  that  the  test  substance  forms  a 
layer  in  the  tube  between  z  =  0  and  z  =  L.  (Figure  1)  Above  the  test  substance 
we  have  a  brass  disk  of  thickness  U  and  below  the  test  substance  we  have  a  water- 
filled  section  S  extending  from  z  =  -  h  to  z  3  0.  Let  o0  -  density  of  the  water, 
cQ  =  velocity  of  sound  in  the  tube  between  z  *  -•  h  and  b  =  0,  Z0  3  specific 
acoustic  impedance  of  the  water-filled  section  of  the  tube.  P+°  and  P_°  are 
defined  as  the  pressure  amplitudes  of  the  incident  and  reflected  waves  in  S, 
respectively.  The  total  pressure  in  F  is  then  given  by: 

# 

135)  P0  =  (P+°‘  e-V5  +  P„°  e+^>z)  €r^et 

where  Te  o  +  3  /co«  The  constanto(0  is  an  attenuation  constant  which  is  now. 
being  introduced.  We  then  have  a  corresponding  equation  for  the  local  velocity 
V0  in  $  so  that  the  specific  acoustic  impedance  Z0  becomes: 


136) 


w 


“  ±  (Jsoo/to) 


'O 


(P+°  e"x°z  +  P_° 

(P+°  e ~Z°z  -  P.°  e"X°z) 


To  determine  the  value  of  Z0  at  z  3  -  h  or  at  the  interface  z  =  0  we  need 
only  replace  Z  in  (136)  by  -  h  or  0,  At  z  =  0  (the  interface)  wc  have 


as  usual. 


In  water,  we  can  neglect  the  attenuation  coefficient so  that 73 0  reduces 
to  jf./c0  as  in  our  former  equations.  Replacing  P_  and  P+.  in  (123,)  by.  P„°  and  P+° 
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we  may  write  P_°/P+°  =  re“3e  where  \i  =  pj  +  jp.R,  r  =  e^P-t  and  9  =  ~2p  as 
before.  Taking  the  value  of  r  in  decibels  we  get 

138)  db  =  20  log  r. 

r 

Numerical  values  of  Z0/o0c  may  be  obtained  by  Smith* s  charts  as  in  the 
foregoing  formulas. 

Now,  let  us  call  oa  the  density  of  the  sample  and  c^the  sound  velocity 
in  the  sample,  P+*  and  P_  1  are  defined  as  the  incident  and  reflected  waves  in 
the  sample  andX',  the  value  of  X-  in  the  sample.  The  pressure  inside  the  sample 
becomes 

139)  «=  (P+*  e~V  A  P_ 1  e4"^iz)  e^£t 

We  have  a  corresponding  formula  for  the  local  velocity  wa  in  the  sample  so 
that  the  specific  acoustic  impedance  Z1  of  the  sample  becomes 


liiO) 


(P+1  e“V  +  P,*-  e^*2) 
(P+1  e~TiZ  -  P_*  e+T«z) 


where  Z  varies  from  2  =  0  to  t  «  L» 
obtained  by  setting  z  ®  0  in  (lliO). 


The  impedance  Z10  at  the  interface  is 
We  have 


The  impedance  at  the  brass  washer  may  be  obtained  by  setting  z 
(lliO),  We  have  for 


=  L  in 


1*  -  e^1 

,  P+1 

lh2) 

ziL  D  (oec^/cx  )  - - - 

1- 

P  + 

Solving 

(lii2)  for  P-  /?+*  we  have 

ih3) 

r-’  .  0-Z\lX.hZ%  ■  3e°l 

P+1  je  Oi  +  \  Z1L 
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From  (139)  we  see  that  the 


lhli) 


F 

lL 


(P+1  e“riL  +  P..1 


and  the  corresponding  velocity 
(118)  to  (139)  so  that 


116) 


XL 


1  dTj 
jec^  dz 


pressure  at  z  ~ 

e*V)  e*** 

at  2  =  L,  w  ^  is 


L,  P^  is  given  by 


obtained  by  a] plication  of 


=  hi.  (P+i  e^i1  -  e^i1)  e^et 

Oec^ 

The  brass  interface  may  be  regarded  as  a  perfect  reflector  so  that  at 


lh6)  v1L  -  ° 

Thus,  from  (Hi5) 

lh7)  P+x  e"\L  -  P_a  e^xL  =  0  or 


e-2H&j  L 


Returning  to  equation  (llal)  we  may  write 


lhR) 


(JfiOiAi) 


1+  e~^il. 
1-  e  -^iL 


Dividing  both  sides  of  (llid)  byec^L,  where  L  is  the  length  of  the  sample,  and 
multiplying  numerator  and  denominator  by  we  obtain 

11.9)  Zl%o,L  -  — 

V 

which  reduces  tc 

Z-in  * 

150)  /(:o,L  *  -i-  coth  X/  L 

1  ZtL  1 


*  e'V 

e®i  E  _  e  L 
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On  setting  Z^q  c  R^q  +  andt/-^  ®  +jp£  and  substituting  coth  Jz  <=  -j  cot  z 
we  have 


1^1)  Rio/ea^ 

2  2 

(J-2L  tanh  p^L  (1+  tan^  p?L)  +  p^L  tan  P2L  (1-  tanh  p],L) 

82  -  — -  - . -  -r-  1  -  -  ■■■  — . —  -  ■ 

2  2  2  2  2  2 
Cm-  ^  2^  )  (tanh-  ai.L  +  tan  ^L) 


°1L 


p^L  tanh  p^L  (1+  tanc  pgL)  “  \Lo 


U  OcUl  \X~Ij  \±- 


2  2  2  2  2  2 
(m-  -jL  +  |i  2^  )  (tanh”  +  tan  pgL) 


Since  p^  and  p2  are  known,  we  may  calculate  X]n^C]L  from  (l£l). 

Special  charts  are  available  which  plot  the  functions  (l£l),  facilitating 
computations.  Once  (l£l)  are  known,  the  components  of  the  impedance  of  the 
sample  may  be  determined. 
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17.  The  Standing  Wave  Fethod 


The  measurements  of  Z  ur  to  now  have  been  considered  as  performed  by  means 
of  the  pulse  tube  method.  We  had  to  determine  the  complex  ratio 


and  the  physical  quantities  accessible  to  measurements  we re  the  amplitude  r  and 
the  phase  angle  Q  at  certain  interfaces  z  =  0  and  2,  =  T,.  'She  relationship  between 
these  numbers  and  the  numbers  Rpo/koiL  and  XpQ,4  OpL  of  Z,^/ea^L 

cR10^alL  was  then  determined  from  equations  (l5l).or  from  special  charts 

of  equations . (l5l) . 

When  measurements  by  standing  waves  are  performed,  the  fundamental  quantities 
to  be  measured  are  the'  ratio  Pmin/l'max  of  the  minimum  to  maximum  pressure  and  the 
phase  difference  as  obtained  by  measurement  of  the  distance  of  the  pressure  nodes 
from  the  sample  interface.  The  equation  for  pressure  (113)  and  the  relationship 
given  by  (123)  still  applies.  Then  P+e"h-  =  Pn  “  P..  e  K  and  (113)  can  be  rewritten 


19) 


P  =  po  (e  u-jez/c  +  e-(p-jez/c)';  e±jet 


P0  is  a  complex  quantity  whereas  the  pressure  measured  at  different  points 
in  the  tube  is  real.  To  establish  a  relationship  between  these  two  quantities, 
we  must  deal  with  the  absolute  value  of  P, 'which  is  a  real  number  defined  by 

|  p|  -  y/??* 

where  P#  is  the  complex  conjugate  of  P,  Then,  from  (152) 

153)  |  P  |  -  P0  I  U-jWc  +  e-(n-o«/c)| 


Letting  p  =  pp  +  jp2  and  noting  that  2  cosh  x  =  ex  +  e"x  we  can  write  (153) 
in  the  form 

l5h)  |P  |  =  |  2PC|  | cosh  (pp  -j  (ez/c  -  p2))| 

~  |  ^  r0  j  |  COS  «/c  P2 ) 

-  j  sinh  pp  sin  (e  z/'c  -  ^2 )  | 

Now,  the  absolute  va]ue  of  P  is 

155)  |P  \  T  1 2P0  1/ cosh2  pp  cos2  (ez/c  -  p2)  +  sinh2  41  p  sin2  (ez/c  -  p2) 
which  can  be  written  (cosh2  pp  =  sinh2  pp  +  l') 


156) 


I  p  i  H2  |\/sinh2  pp  +  cos2  (sz/c  -  p2) 
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We  see  that  j  P  J  is  a  minimum  if  cos2  (ez/c  -  p2)  °  0  which  is  satisfied  by 

157)  P2  -e  z/c  «  (n  +  $)*,  n  =  °,  1,  3  .  .  .  . 

The  maximum  of  | P  |  occurs  if  cos2  (ez/c  +  p2)  “  1  which  corresponds  to 

158)  p2  -  ez/c  s  ri^1,  n=  0,  1,  2,  3  »  *  * 

Thus  we  have 

159)  (P  (maximum  **  2|P0jcosh  p-j. 

j  P  (minimum  =  2jp0jsinh  p-^ 

For  a  perfect  reflector,  pi  «  0,  j  P  (minimum  °  0,  j  P  (maximum  *  2 j  P0  j. 

In  other  cases,  the  ratio  * 

I  p  ”inj/ |p  max  | 

is  experimentally  known  and  we  have 


160) 


P  min  |y/j  p  mflx  j  _  tanh  ^  or 

-1  j  P  min  \j  J  p  max  J  or 


if 


p^  =  tanh 
P  min  j  j j 


P  max  j  =  (3  we  nave 
p^  =  tanh  “i  p.  and  thus 


161)  px  *=  (3  +  p  V3  +  P  */5  +  .  .  .  +  (3  n/n  +  .  . 

2 

which  gives  pj,  by  successive  convergent  approximations  for  (3  <L,  a  condition 
always  satisfied. 


The  problem  can  also  be  treated  without  using  complex  numbers.  Consider 
the  local  displacement  caused  by  the  acoustic  waves  in  the  tube.  It  is  given  by 

16?)  h  =  a  cos  k  (ct  -  z)  -  r  cos  k  (ct  +  z) 

=  (a  +  r)  sin  k  ct  sin  k  Z  +  (a-r)  cos  k  ct  cos  kz 

where  a  is  the  amplitude  of  the  incident  wave,  r  the  amplitude  of  the  reflected 

wave.  The  motion  can  be  regarded  as  due  to  the  superposition  of  two  stationary 

waves  of  amplitudes  (a  +  r)  and  (a-r)  and,  according  to  the  definition  of  the 
pressure  P  in  terms  of  the  displacement  h 


K 
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4 


163) 


«  o  c 
0 


dhi 

3t 


c 

0t 


we  have 

l6li)  P  «  P^  +  Pr  =  a0kc^  |^(a  +  r)  cos  kct  sin  kz 

-  (a  -  r)  sin  kct  cos  kz^J 


where  k  =  2Hv/c,V“  being  the  frequency  of  the  sound  velocity  c,o0  the  average 
water  density.  We  have 


165) 


P  min  a  -  r  1 

■  ■  .  -  ■■  «  — — ^  —  m  -  -  -  - - - - -  . 

P  max  a  +  r  Standing  Wave  Ratio 

The  reflection  coefficient  r/a  of  an  imperfect  reflector  is  expressed  as 


166) 


r_  d  SWR  -  1 
a  °  SWR  +  1 


For  example,  if  SWR  =  2,  -~ 


1 

3 
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18. 


Illustrative  Example 


11,  12 


Consider  a  set  of  measurements  to  be  made  by  the  pulse  tube  method.  The 
freq:  „ncy  range  is  from  2,000  cps  to  20,000  cps  and  there  are  a  minimum  of  10 
waves  per  pulse  trajji.  The  characteristic  impedance  of  fresh  water  is 
00c0  =  1.5x105  gm/cm-  sec:  the  thickness  of  the  sample  is  1.055  cm  and  its 
density  is  op  =  1.18  gm/cm3.  The  operating  frequency  is  15,000  cps.  We  would 
like  to  determine  the  characteristic  impedance  of  the  sample,  its  attenuation 
coefficient  and  the  velocity  of  sound  in  the  sample.  Measurements  indicate  that 
the  reflection  coefficient  is  0.91:5  and  the  phase  change  -  31«1°. 

At  a  minimum  frequency  of  2000  cps  the  wavelength  per  train  of  ten  waves  is 


X  =  —  n  =  |222  x  10  =  25  feet 
V-  2000 

The  minimum  theoretical  length  of  t'he  tube  is  therefore 

L  =  —  =  12.5  feet 
2 

Choice  of  this  length  implies  that  there  is  no  time  delay  between  the  instant 
the  wave  train  leaves  the  oscillator  and  the  instant  it  returns.  Under  a;tual 
conditions  there  is  a  finite  time  lapse  between  projected  and  reflected  wave 
trains,  and  the  tube  is  generally  made  longer  than  theoretical  length  to 
accommodate  this  time  difference.  We  shall  choose  a  tube  factor  F  =  0.75  so  that 

L  »  FA  =  0.75  x  25  =  19  feet 


is  the  actual  length  of  the  tube. 


If  we  choose  our  axes  at  the  sample  interface  we  have,  from  (122) 

1  Iz 

Zoo  _  r+ 


°oco 


1  - 


where 


L  oo 
a_c_ 

KJ  U 


=1/00=u  4 


i  - 


p_ 

jv  and  --  =  x  +  jy  ®  \  . 


p 

The  complex  ratio  ~  is  termed  the  reflection  coefficient  and  consists  of 

*  + 

an  amplitude  and  a  phase  change,  which,  for  this  problem  are  0,91:5  and  -  31.1° 
respectively.  Hence,  we  may  rewrite  (122)  as 

^  =  Lld 

^  00  1  -  T, 

where -q  *  x  +  j  y  1  0.91:5  e~3  (31*1  )t  Thus, 
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x2  +  y2  «  (0.9b5)2  -  0.893 
167)  and 

tan  (-  31.1°)  =  •%-  *  -  O.60I4 

X 

solving  equations  (167)  we  have  x  °  .808  and  y  =  -  .L188  so  that  -q  =  0.808  -  j  O.L188. 
From  (130)  we  have 


2  2 
1-y  -x 

( 1 —x ) 2  +  y2 


.397 


+  jv  =  0.397  -  j 


v 


?y  b 

(1-x)2  +  y2 


-3.55 


3.55. 


The  same  result  can  be  quickly  obtained  from  a  Smith  Chart  (Figure  3). 
Equations  (131)  represent  circles  in  the  (x,  y)  plane  with  centers  at  (,28b,  0) 
and  (1,  -  .282)  respectively  with  radii  of  0.?l6  and  0-282,  To  determine  u  and 
v  from  a  Smith  Chart,  draw  a  vector  of  length  r  =  0.9h5  at  an  angle  of  -  31.1° 
with  the  x  axis.  The  head  of  the  vector  r  lies  at  the  point  of  intersection  of 
the  two  circles  u  and  v.  Values  of  p.  =  O.liO  and  v  *  -3.50  may  be  read  directly 
from  the  chart  without  further  calculation.  'Thus 


oo'jf*  0-397  -3  3.55  - 

o  0  0  0 


Since 


we  have 


J10 


%L  =  ?io_  ::  ^ oc( 

°leL  Vo  °l£l 


1.5  X  10- 


o-j^eL  1.18  x  2  x  l5oCO  x 


T7o?i  (°-397  - j  3-ss) 


The  function 


Z-)  n  Z 


—  =  42-  X  J.285  «-0.5l— 


o1eL  a0cQ 


(151)  |2.  -  o.5l  -  3  h.56  - 12_  ♦  3  |2- 
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is  a  function  of  the  parameters  an^  U2L  anc*  ls  avallallle  ln  chart  form 
(U.S.  Navy  Underwater  Sound  Laboratory  Chart).  A  typical  chart  is  shown  in 
Figure  h.  Entering  the  chart  with  values  of 

=  0.5l  and  =  -  hr56 

o-jeL  o-j_cL 

we  find  that  PqJ,  =  .025  and  n2L  =  .1.55. 

Values  of  p^L  and  p2L  obtained  from  the  NAVU/JTRSOVNDLAB  chart  may  be 
verified  by  substituting  ji^L  =  0.025  and  p2L  =  0.h$$  in  (l5l)  and  calculating 
Rio/oieL  and  X^q/o^eL.  Thus 


RlO/b-^L  c 

A6$  tanh  .025  (1*  tan2  .1.55)  +  ,025  tan  ,1.55  (1-  tanh2  .025) 
(.0252  +  .1.552)  (tanh2  .025  +  tan2  ,Ii55) 

Rio/aieL  =  0.528 

*1  = 

.025  tanh  .025  (1+  tan2  .155)  -  .ifeo  tan  .1.55  (L-  tanh2  .  .025) 
(.0252  +  .h552)  (tanh2  .025  +  tan2  ,h$$) 


Xio/oieL  =  -  lt.l.li9 

R-i  q 

These  values  agree  within  lx%  of  the  chart  values  ofa~^ 


°.5l  and^ 


-  -a .56. 


Thus 


PXL  -  0.025;  ji1  =  =  .0237 


M-2L  -  0.1.55;  P2  =  =  .li31  cm"1 


and,  since 


Y*  zio  <y  .  6 
®  10  =a—T  1  +  j  CT  =  ^1  +  5  ^2 
11  1 

we  have  | 

oC.  i  =  Pi  =  .0237  =  attenuation  coefficient, 


— -  e  ^2  3  »^31  so  that  the  velocity  of  sound  c-^  in  the  material  becomes 
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2*x  1$000  ,l0  , 

- - -  =  2.18  x  1CK  cm/sfcc 

.1j31 


and  the  characteristic  impedance  of  the  sample  is  therefore 


OnOj  »  1.18  x  2.18  x  10*  =  2.3'?  x  10$ 
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PULSE  TUBE  SCHEMATIC  FOR  THE  MEASUREMENT  OF  THE  SPECIFIC  IMPEDANCE 
OF  A  BRASS  DISC  IN  A  WATER  FILLED  TUBE 


FIGURE  2 


SMITH  CHART  REPRESENTATION  OP  THE  FUNCTION  /=  (l+£)/(l-£) 
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